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ANNOTATION 


The present book represents an attempt at uniform consider- 
ation of the various sines (the circular and hyperbolic sines, 
the lemniscate sine, and the sine amplitude of Jacobi) as special 
cases of the so-called generalized sine—all of them functions 
defined as the inverse of some integral. It requires a certain 
amount of mathematical training and is intended for sufficiently 
prepared readers who have mastered mathematical analysis as 
taught in university courses. 
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PREFACE 


Of the functions studied in elementary mathematics, the 
trigonometric functions stand out by virtue of their geometrical 
definition. Without stopping for inconsequential variations of the 
same idea, we may say that the sine and cosine are introduced 
as coordinates of a point on the unit circle. The independent 
variable is treated as an angle or arc of a circle. 

In the present book, we shall show how it is possible, by 
beginning with other curves (such as the equilateral hyperbola 
or Bernoulli’s lemniscate (a curve having the form of a figure- 
eight), to define interesting and important functions analogous 
to the trigonometric functions, similar to theminsome respects 
but possessing certain new characteristics. These functions are 
called respectively hyperbolic and lemniscate functions. In 
analogy with them, we shall refer to the trigonometric functions 
as ciyculay functions. 

We shall consider all these functions as special cases of the 
generalized sine, that is, the function inverse to the function 
defined by an integral of the form 


y 

{ dz 
Xe ao 
: V 1+ mz? + nz! 
0 


Here, the circular sine corresponds to the casem=— Il, n=O; 
the hyperbolic sine corresponds to the case m=—1,n=-0; the 
lemniscate sine corresponds to the case m=0, n=-—\, If 
m= —1—k? and nk? (for 0<& < 1), we obtain for the function 
inverse to that defined by the integral the so-called sine ampli- 
tude of Jacobi. The problem of the pendulum leads to a study of 
this last function, 


x PREFACE 


To study the properties of all these functions in a uniform 
way, we need first of all to define them as functions of a com- 
plex variable and then establish an addition theorem for them. 
This path was already taken at the end of the eighteenth century 
by the young Gauss in his Tagebuch. 

The reader is assumed to have a familiarity with the ele- 
ments of analytic geometry anddifferential and integral calculus. 
The necessary material on integration in the complex plane will 
be given in the present book though proofs will be omitted. 

The ultimate purpose of the book is to acquaint the reader 
not possessing an extensive knowledge of the theory of functions 
of a complex variable with the simplest representatives of the 
class of elliptic functions, namely, lemniscate functions and the 
somewhat more general Jacobian elliptic functions. 

In conclusion, we warn the reader that this book is not in- 
tended for light reading. He must read it with his pencil in his 
hand, 


The Author 


CHAPTER 1 


Geometric Definition of Circular, Hyperbolic, 


and Lemniscate Functions 


1. The circular functions are defined and their properties 
are discussed in high-school mathematics classes. Here, we 
shall briefly expound familiar facts in a form that will be con- 
venient for us to bring out the similarities and dissimilarities 
sharply. 

Let us begin with the unit circle 


xp yh L (1) 


The arc length from the point C toa point A (see Fig, 1) is 
defined up to an integral multiple of the total circumference of 
the circle, namely, 2nk, Here, k is the number of complete 
revolutions about the coordinate origin described by a point on 


FIG. 1. 
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the circle as it moves from C to A. We take the counter- 
clockwise direction as positive. 

The coordinates x and y of the point A are independent of the 
number of complete revolutions made. Therefore, they are 
periodic functions of ¢ with period 2n. By definition, 


x= cost, yo sine: (2) 
From this it follows that, for arbitrary ¢, 
cos?#-+sin?t= 1. (3) 


It follows from the symmetry of the circle about the x~axis 
that 


cos(— f)==cosf, sin(— ¢)= — siné. (4) 


If we also use the symmetry about the bisector of the first and 
third quadrants (by virtue of which the abscissa and ordinate of 
any point on the unit circle are interchanged), we see that 


COs (> —~ t) == Sinf, sin ($—<) = cos f. (5) 


Each of the formulas (5) enables us to obtain a definition of 
either of these two functions in terms of the other; for example, 
we can define cos ¢ in terms of sin ¢ (assumed already defined). 

As consequences of formulas (4) and (5), we can derive all the 
simple ‘treduction’? formulas. Thus, for example, we obtain the 
following pairs of formulas; the derivation of each pair after 
the first uses the preceding formulas: 


nN 


cos (=; ae t} = COS E 


sin( + t) = sin 


( | —sin(— d= — sint, 


Ja 


a | == cos(— f)= cos f; 


2 
costa — f)== cos |= of 5 —O\|=— sin (= —ft\)=-— cost 
. _ 12 } aaa pi ’ 
Sin = t= sin] 5 +(5— \|= cos (F — thas sing 
7 sé 
cos (-- f)== cos[1 — (— fF) = — cos(-— t)= -- cos f, 
sin (a+ f) == sin{[a-- (— #] = sin(— f) = — sin; 
cos (2m -+- f) = cos [1 4-(1 + 1)) == — cos(a1+-f)=—cost, 


sin( 22 4- ¢) == sin[a+-(1-+- 4) == — sin(a 4 f) =sinl. 
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We note that the second pair of formulas can be derived imme- 
diately from the symmetry of the circle about the y-axis and 
that the third can be derived from the symmetry about the co- 
ordinate origin. The fourth pair expresses the periodicity, 
mentioned above, of the circular functions. 

Defintion (2) and equation (1) enable us to answer the question 
regarding the number and nature of the solution of equations of 
the form 


sinf=-a@ or cost=—a, 


where « is a real number, 
For example, consider the equation 


sint =a. (6) 
We may assert that, for |a| > 1, this equation has no solution 


at all. For |u| <1, there exists at least one solution in the 
interval 0 <¢< 2a. Specifically, 


if ==, the only solution is t=; 


if 0<a< 1, there are two solutions ¢, and 4,; 


furthermore, 0 <t,< 5 and > <ty=n—t <1 
if oa-0, there are two solutions, namely, ¢, 0 and t, ==; 
if —Il<u< 0, there are two solutions ¢, and 4); 
furthermore, 1<t, < aa and s < fy =e bn — ty, < 20; 


3m 


if a= — 1, there is only the solution t=-. 


Of course, in the cases of a equal to 1 or —1, we can speak 
of two solutions that coincide. If ¢,#4¢,, the corresponding 
values of the cosine are negatives of each other: 


cost, = — cos fly. 
Geometrically, solution of equation (6) reduces to finding points 


on the circle (1) of given ordinate y =« and the evaluation of the 
lengths of the arcs terminating at the points found (see Fig. 2). 
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Allthat has been said in the present section is well illustrated 
on the graphs of sintandcost, The construction of these graphs 
depends somewhat on the definitions accepted (see Fig. 3). 


FIG. 2. FIG, 3. 


We note that we might interpret the parameter ¢ in the defini- 
tion (2) geometrically not only as the arc length CA but also as 
twice the area of the sector OCA: 


t=2 area OCA. (7) 


Here, we need to regard the sector as produced by a moving 
radius displaced from the position OC to the position OA. If the 
radius undergoes k complete rotations in the process (where k 
can be either negative or positive (according to the direction of 
the rotation), the area swept by it will be equal to ka. Thus, we 
have the term 2kn in the definition of ¢ given by formula (7). 

2. Now, instead of the unit circle, let us consider the unit 
equilateral hyperbola 


xe— y= l. (8) 


Equations (1) and (8) resemble each other much more than do 
the curves representing them (cf. Figs. 1 and 4). In fact, this 
is putting it mildly. Equations (1) and (8) differ from each other 
Only by the sign in front of y?.. On the other hand, the curves 
shown in Figures 1 and 4 are not the least bit like each other. 
Yet, it is the equations that we must trust rather than our eyes. 
The similarity in the equations entails a far-reaching similarity 
in the properties of the curves, which we discover when we 
compare the circular with the hyperbolic functions. 
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FIG. 4. 


Let us denote by ¢ twice the area of the sector OCA: 
t=2 area OCA, 


We shall consider the area of OCA positive if the rotation from 
Oc to OA is carried out in the counterclockwise direction and 
we Shall consider it negative in the opposite case. If the point A 
describes the right branch of the hyperbola shown in Fig. 4 
as y increases from —oo to + co, then ¢ increases also from 
—« to +oco, The coordinates x and y of the point A can be 
regarded as single-valued functions of f. We define 


x=cosht, y=sinht. (9) 
Here, the abbreviations ‘‘cosh’’ and ‘‘sinh’’ are for the expres- 
sions consinus hyperbolicus (hyperbolic cosine) and stmus hyper- 
boltcus (hyperbolic sine). 
It follows from the definition that, for arbitrary +, 


cosh7¢—sinh7t= 1. (10) 


Using the symmetry of a hyperbola about the x-axis, we 
conclude that 


cosh (—t)=cosht, .. sinh (~t)=—sinht (11) 


(cf. the graphs of cosh ¢ and sinh ‘in Fig, 5). 
The equation 


sinh (t) =a (12) 
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has one and only one root for arbitrary real a. The equation 


cosh t=a (13) 
has no root at all ifa< 1, one root (namely, ¢=0) ifa—1!, and 
two roots ¢, and ¢, (-==—4t,) ifa@>1. 


Geometrically, the solution of equations (12) and (13) can be 
reduced to finding points of given ordinate y» or given abscissa 
x > !on the right branch of the hyperbola (10) and to the evalua- 
tion of the areas of the sectors corresponding to the points 
found (Fig. 6). 


J 


FIG. 5. FIG. 6. 


3. Before defining the lemniscate functions, we note that it 
is possible to givetothe circular functions a geometric interpre- 
tation different from the interpretation given in section l. 
Specifically, let us look at the circle of unit diameter tangent 
to the x-axis at the origin and lying in the upper half-plane: 
x?+ y?— y= 0 (see Fig. 7). Obviously, the chord OA connecting 


FIG. 7. 


the end points of the arc O4 of length ¢ is equal to sin ¢. There- 
fore, it would be possible to construct the theory of circular 
functions by defining 
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sint-=OA, t-= length OA: = 4\woshe DR (14) 


It follows from this definition that sin ¢ increases from 0 to 1 
as ¢ increases from 0 to «/2. As t increases beyond 1/2, the 
Sine decreases and vanishes as tn. (Under the conditions 
stated, nm is the circumference of the circle.) 

Let us complete our definition by agreeing that sin ¢ changes 
Sign every time that it passes throughthe value 0 as ¢ increases. 
Then, we need to take the length of OA in formula (14) with the 
minus sign for n</¢< 2n and with the plus sign for2a <f < 3n, 
etc. Only under such a supplementary convention will our new 
definition coincide with the preceding one, and only then, in 
particular, will the period of the sine be 2n (rather than x as 
would be the case if we considered the length of the chord OA 
positive for all possible values of the arc OA, which would then 
differ by integral multiples of x). 

The second circular function, the cosine, is defined in terms 
of the sine by the equation 


cos f= sin (5 — t). (15) 


4. Let us look at Bernoulli’s lemniscate shown in Fig, 8. 
This is the curve consisting of all points in the plane such that 
the product of their distances from two given points, called the 
foci of the lemniscate, is constant. In the present case, the 
points (1/2, 1/2) and (—1/2, —1/2) are chosen as the foci and 
the constant distance product is 1/2. Therefore, the equation 
of the lemniscate is of the form 


nn fe 


\ 


i 
%| 


or, after squaring both sides and simplifying, 
(x? 4 yt)? — Dey = 0. (16) 
Thus, Bernoulli’s lemniscate is a fourth-order algebraic curve. 
Let us introduce polar coordinates r and y, taking the x-axis 
as the polar axis. Then, 


X= rcosfi, Yue rsing. 


Thus, the equation of the lemniscate in polar coordinates 
becomes 
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r? = sin 29. (17) 


The variable , attains its maximum value, namely,1 at p= 1/4. 

It follows from equation (17) (or (16)) that the lemniscate 
is symmetric about the coordinate origin. Furthermore, it is 
symmetric about each of the diagonal lines bisecting the angles 
between the coordinate axes: if x =aand y=»+ satisfy equation 
(16), so do the pairs (— a, —bd), (6, a), and(—b, —a). 


FIG. 8. 


5. Suppose that a movable point ™ starts at the origin and 
moves in such a way astodescribe that portion of the lemniscate 
lying in the first quadrant, moving in the positive (i.e., counter- 
clockwise) direction. Since it is being displaced along an arc 
tangent to the y-axis when it first returns to the origin, it is 
natural to require that it be moving along an arc tangent to that 
axis as it enters the next (third) quadrant in the course of its 
motion. Accordingly, we must take the clockwise direction for 
the positive direction of motion of the point as it moves around 
that portion of the lemniscate in the third quadrant. When the 
point is finishing its first complete trip over the entire lemnis- 
cate and is beginning a new one, the point leaves the third and 
enters the first quadrant along an arc tangent to the x-axis. The 
motion thus described is then repeated. 

Let us denote the length of the entire lemniscate by 20. Thus, 
the length of the arc OMC is equal to w/2, Just as in the case 
of a circle, we shall assume that the variable point M moves 
from an initial point O of the arc to its terminal point A after 
some number of complete courses over .the entire lemniscate 
in either the positive or the negative direction. Thus, the length 
f of an arbitrary arc OA is determined only up to an integral 
multiple of 2. Therefore, r, the length of the chord OA, isa 
periodic function of ¢ with period 2, 
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We complete the definition of this function by agreeing that it 
changes sign every time it passes through the origin with con- 
tinuous increase in f¢, 

In analogy with the circular sine, we shall refer to the func- 
tion that we have just defined for all real values of ¢ as the 
lemniscate sine and we shall denote it as follows: 


OA==r-—=slt. (18) 


Here, the letters ‘‘sl?’ are the initials of the expression sinus 
lemniscaticus. We note in particular that 


si 5 =. (18’) 


From the symmetry of the lemniscate about the origin and 
the convention adopted regarding sign, it follows that sl ¢ is an 
odd function (see Fig. 8): 

sl(—- f) = — sit. (19) 
From the symmetry about the bisector of the first and third 
quadrants, we have (see Fig. 9) 


sl(w — f)= slit. (20) 


A fa+t) 


FIG. 9. 
Analogously, from the symmetry about the bisector of the 
second and fourth quadrants, we have (see Fig. 9) 
io t)==— sll. (21) 
Actually, (21) is a simple consequence of equations (19) and (20): 


si(m + 4) = sl[w —(— A) = sl(--H#=>— slt. 
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The lemniscate cosine is defined in terms of the sine by 


clt=sl (5 — t). (22) 


\ 


From this, it follows in particular that it, too, is a periodic 
function with period 20, However, in contrast with sl ¢, the 
function cl ¢ is an even function. To see this, note that, by using 
successively (22), (21), (19), and again (22), we have 


w 


el(—#) = s1($ sere =—sl(—3 +4] eae ($—t)= = cle. (23) 


The reader can easily verify the following identities: 
ci(m--— ft) = —clf, cl(w + O-= —clt. (24) 


The definition of the lemniscate functions and equation (16) 
enable us to answer the question as to the number and nature 
of solutions of equations of the form slf—a or cl t==a, where 
ois a real number. 

Here, we have a perfect analogy with the case of circular 
functions (cf. section 1). 

For example, let us consider the equation 


sl f=. (25) 


We may assertthat, for ]a| > 1, it has no solution at all. For 
|a|<< 1, there exists at least one solution inthe intervalO<t< 2. 
Specifically, 


if a= 1, there is the single solution ¢= 33 


if U<a< 1, there are two solutions, ¢,; and 4; 
these are such that0 <4, < $ and 5 <t=o—t, <0; 
if a=0O, there are two solutions f, = 0 and ¢, =a; 


if —1<a<0, there are two solutions ¢, and 4; 


30 


these are such that wo <f, ee = and < ty = 30 —t, < 20; 


if a == — I, there is the single solution os 
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Obviously, in each of the cases a=1 and a=-— I, we can speak 
of two solutions that coincide. 

We note also that, if ¢, #4,, the corresponding values of cl ¢ 
are negatives of each other: 


cl f; = — cl f. 


Geometrically, solution of equation (25) reduces to finding the 
points of intersection of the lemniscate (16) and the circle 
x?+ y*--q?, to the choice of those points on the intersection for 
which r=-OA has the same sign as a, and, finally, to the cal- 
culation of the corresponding arc lengths. 

All that has been said in the present section is illustrated in 
the graphs of sl ¢ and cl ¢, which, in appearance, resemble the 
sine and cosine graphs (see Fig. 10). 


CHAPTER 2 


Generalized Sines 


6. Let us now show that a uniform approach is possible to 
the study of the sines (and with them the cosines) defined 
geometrically in Chapter 1. With this end, we shall consider 
each type of sine as a function inverse to some integral. 

Let us begin with the circular sine, At the end of section 1, 
it was noted that the variable t can be interpreted as twice the 
area of the circular sector OAC (see Fig. 1). If x and y =sin? 
are the coordinates of a point A, then x= |) “T— yw? and, con- 
sequently, 


area A OAK = sy VI. 


The area of the curvilinear trapezoid VUCAK is calculated from 
the formula 


y y 
area OCAK = f dy== | V1 — Fay, 
6 0 


Integrating by parts, we obtain 


‘i \l — wdqn= VY 1— y? y +. / Jie = 


y 


y 
S34 oa Pe a 


0 


14 THE REMARKABLE SINE FUNCTIONS 
from which we get 


y y 
§ _ ees I 
| Vi-van=sVi- yt f- dn 
0 


. 2 Vt—y : 
Consequently, 
Z 
S= area OCA=area OCAK—area * OAK = a | meee ; 
~ 28 Vi—7 
or 
: ea 
= 25 — aa ae (26) 
; / VI? 


0 


As the variable y increases from —1 to 1, the variable ¢ in- 
creases from —n/2 to 1/2. We can calculate 1/2 from the 
formula 


ay I 
== — = 1.5708. 27 
| Vi—7Y? 2 ( ) 


Conversely, as ¢ increases from —n/2 to 2/2, the variable 
v==sinf increases from —1to1. Thus, we can regard the func- 
tion y-=sin¢ in the interval (—n/2, 1/2) as the function inverse 
to the integral (26). As before, the function cos ¢ is defined by 
the formula 


na 
¢=sin(>—2). 
cos sin (5 


7. Let us now look at the hyperbolic sine sinh ¢, In defining 
it geometrically (cf. section 2), we interpreted the variable ¢ 
as twice the area of the sector OCA (see Fig. 4). We have 

S=area OCA=area OCAK — area OAK. 


By using the equation for a hyperbola x? — y?=1, we find 


area OAK an0 xy=s Vit+y?y 
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and 
y y 
area OCAK = [ean= [ ViFFan, 
0 0 


Integration by parts yields 


l 


s y 
° ——s = ——; = a n? 


6 


so that 
y y 
aera laos Lf an 
area OCAK = l+rndyn=——V1+y*y+— —— 
ny > | ad Vi+ 7% 
Consequently, 
J, 
Pe. Ae dn 
S=— 
2 V 14-7? 
or 


y 
dy 
f Vitn 


As the variable y increases from — cc to + .w, the variable ¢ 
also increases from —co to + oo since 


[ ay rae 
4 Vi + 7? 


Conversely, as f increases from —cc to -+ oo, the variable y 
also increases from —cc to + oo. Thus, we can regard the 
function y = sinh¢ on the entire real axis (—oo, +o) as the 
function inverse to the integral (28). 
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8. In what follows, we shall choose a procedure of exposition 
in which the integrals (26) and (28) willbe the original functions, 
and the functions sin/ and sinh? will be secondary, that is, 
obtained from the preceding. We shall define the latter functions 
as the inverses of the integrals. Of course, when such defini- 
tions are made, the integrals themselves can be regarded as the 
functions inverse to sin ¢ and sinh? and we can introduce the 
following notations for them: 


v 


==aresiny and ¢ = | aan = ==arsinh y. (24 
é +4 


y 
t= f — 
yee 


In the first case, the syllable ‘‘arc’’ is appropriate because 
twice the area of the sector OAC in the unit circle is measured 
by the same number as is the arc AC, In the second case, the 
syllable ‘tar’’ stands for ‘‘area’’. In the case of the hyperbola, 
however, there is no such simple relationship between the area 
of the sector OAC and the arc AC as in the case of a circle. 

Finally, formulas (29) do not as yet provide us with new 
means for calculating the corresponding integrals. For example, 
the first of them simply asserts that the integral is the function 
inverse to sin f, but we have decided to define sin ¢ itself as the 
function inverse to the integral! 

However, for the second of the integrals (2), there is an 
analytic representation that does not depend on our definitions. 
Specifically, it is given by the following formula, the validity 
of which can be shown by differentiating: 


y 


: dy Soe mEEET 
fas | View = = In (y “t Vi le y’). (30) 


This enables us to obtain a formula for the nite) Matsa 
y = sinh ¢ without needing to resort to the integral (38) every 
time we need to speak of the hyperbolic sine. Specifically, it 
follows from (30) that 


so that 
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If we subtract the second of these equations from the first, we 
obtain 


ese" 


y=sinht= (31) 


Thus, it would have been possible at the very beginning to 
define the hyperbolic sine from formula (31) as half the differ- 
ence of two exponential functions. 

Then, formula(10)for cosh ¢ yields 


t -\2 rt ~f \2 
~h2 e—eé ete 
cosh? ¢ = | +(£43— = pon : 


Since (ef+e7')/2 >0 and cosh ¢ is also positive (we recall that 
x =coshf is the abscissa of a point on the right branch of the 
hyperbola (Fig. 4)), we have 


-t 
cosht ag (32) 


This formula defines the hyperbolic cosine as the average of 
the exponential functions e! ande’‘. In Chapter 5, we shall see 
that the circular functions can also be expressed in terms of 
exponential functions though only with imaginary arguments of 
the exponentials. 

9. Let us turn to the lemniscate sinesl] ¢, where the variable 
/ is treated as the arc lengthofthe lemniscate, Let us use polar 
coordinates for the integral expression for ¢, In polar coordin- 
ates, the differential of the arc length is given (see Fig. 11) by: 


dt= Vir dpe + dr = Vac (ty 42 Pare 
vy 


In the present case, equation (17) yields 
C= 5 arcsin (r*). 
From (29) (with yr’), we obtain 


Ek Ane ees Sete 
dr Vl=r 
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Consequently, 
r’ ar 
‘= + ldr= 
Vo ay : Vi-r 
so that 
dt 
=—/ ==. 33 
t | Vt (33) 
As r increases from —! to -++1, the integral (33) increases from 
—w/2to o/2, where »/2 has the value 
] 
= | ve = 1.3111. (34) 
/ — ri 


Conversely, as ¢ increases from — w/2 to o/2, the variable r 
increases from —1l to 1. Thus, we can define the function 
y=sl¢in the interval (— o/2, »/2) as the function inverse to the 
integral (33). As before, we shall define the function cl ¢ by the 
formula 


clf=si(5—¢). 


FIG. VN. 


10. Let us compare the definitions of the three sines given 
in this chapter: (1) y=—sinf (the circular sine), the function 
inverse to the function 


y 


lve 
l—1 
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and defined in the interval (—x/2, 1/2), where 


1 
xf an 
20 J Vi=q’ 


(2) y = sinh ¢ (the hyperbolic sine), the function inverse to 


y 
lay 


0 


and defined on the entire real axis; (3) y = sI ¢ (the lemniscate 
sine) the function inverse to 


y 
= 
° —n’ 
- Vin 


and defined in the interval (—w/2, 0/2), where 


I 
Oo _ | ay _ 
2 Vi-W 


e 
0 


Obviously, all these cases will be included if we succeed in 
studying the function y==s(f) inverse to an integral of the form 


y 


ey dy 35 
‘= | ——— 
| V 1+ mn? + nm! i) 


where m and 7 are any real numbers, Let us agree to call such 
a function a ‘‘generalized sine.”’ 

For m=—1 and 2=0, we have the circular sine; for m= | 
and n=0, we have the hyperbolic sine; for m=0 andna=. — 1, we 
have the lemniscate sine. 

If m=—(1+ a?) and n=k?, where 0<4< 1, the polynomial 
1+ mr’?-+ ny takes the form (i — n’)(1 — 4%). In this case, the 
function s(f) is.the Jacobian elliptic function with modulus a, It 
is called the sine amplitude and is donated by sn(f, &) or, more 
briefly, sn ¢. 

The function sn ¢ is encountered, for example, in the prob- 
lem of the pendulum. This problem consists in study of the 
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oscillations, in a vertical plane, of a heavy sphere of mass m 
at the end of a fine thread of length / (see Fig. 12), Suppose 
that this thread with the attached sphere is held at an angle 
& from its equilibrium position and that the sphere is then 
released with zero initial velocity. 


FIG. 12. 


When the string is passing through the position that makes 
an angle § with the vertical, the velocity of the sphere is 
v ==! di/at and its kinetic energy is 


Yy at 


2 [2 { dQ\2 
mu m (=) 


From the theorem on kinetic energies, this value must be equal 
in the present case to the work done by gravity (we neglect air 
resistance): 


mgB)B = mgl (cos 8 —- cos 9,). 
Thus, 


se (F) = mgl (cos § — cos 6,), 


so that 


Be gee SOE 
2g Vcos8—cos 0, 


= 
| 


(The minus sign is due tothe fact that the angle 0 first decreases 
as ¢ increases). Integrating, we obtain 
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a 7] 0, 
a i d@ a 
22 J V cos 6 — cos 0 V 2 


Let us make the change of variable 


5 () - D, ‘ 
sin 2; ie sin > Ne 


Obviously, n has a simple geometric interpretation, namely, 


We obtain 


) (1 —- 4°) [1 -— sin’ = | 


If we denote the constant 


form 


by ¢, (representing the instant the pendulum passes through 
its equilibrium position), we can rewrite the last formula in the 


V (1 — 9%) (1 — sin 34] 


ae | 

= =e ao Sef ioess § 
¥ rag ae | poe ) 
1 (1 — ny") (| -—— sin: oh ‘\° 


} 


n { dv aa 
ge V cos8—cos 0, 
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and, consequently, 


Thus, we have expressed the quantity n characterizing the dis- 
placement of the pendulum at the instant / from its equilibrium 
position in terms of the Jacobian elliptic function with modulus 
k = sin (u,/2). 

11. It is natural to require that the polynomial 


P(y) = 1 4- mn? + nn 


in formula (35) not be a perfect square. Otherwise, VY P(n) will 
be of the form 1+ q1, and the integral (35) will be expressed 
either in terms of the arc tangent (if ¢ > 0) or in terms of the 
logarithm of a linear fractional function (if g< 0). In each of 
these cases, examination of the inverse function would yield 
nothing new. Therefore, in what follows, we shall assume that 


m—~ 4n =£ 0. (36) 


We also need to distinguish between the following two cases: 

I. The polynomial P()) has no real zeros. This means that 
either the second-degree trinomial 1+ mz-+-nz? has no real 
zeros (that is, m? -4n~<_U) or its zeros are negative real num- 
bers (that is, m’-—4a >> 0, n > 0. m> 0). In this case, P()) >» 0 for 
arbitrary (real) y. Therefore, the integral (35) is defined and 
represents a strictly increasing function on the entire real 
axis. We denote the expression 


| 7 21) 1 
J yl my tn} 


Q 


by 4. Forzn==0, this integral diverges and A-=> For n+ 0, it 
converges and 4 is a finite positive number. Since the, integral 


(35) increases continuously from —A to +4 as y» increases 
continuously from — xo to +20, the inverse function y = s(f) is 
defined in this case in the interval (— A, — A) and it increases 


continuously in that interval from — «.to --.0c, The hyperbolic 
sine (for which 4.-<°o) provides a very simple illustration of 
this case. 


GENERALIZED SINES 23 


II. The polynomial P(n) has real zeros (either two or four 
such). This means that m*—4n >0 and that either 2<0or 
m<O<n. Let au denote the smallest positive zero. Then, —a 
denotes the negative zero of smallest absolute value. Therefore, 
the polynomial P() has no zeros in the interval(—a. a). Hence, 
P(n) does not change sign in that interval; it remains positive 
Since P(U)-=-1>0. [We note that P(y) becomes negative when 
1) (increasing) passes through the value a.] The integral (35) is 
a function that is defined and strictly increasing in the interval 
(—a,a) We denote the expression 


a 


ft} 
J V1 + mr? + ani 


by @ (which is > 0). Since the integral (35) increases continuously 
from —a to a as y increases continuously from —a to a, the 
inverse function y=s(f) is defined in this case in the interval 
(—a, a) and increases continuously in that interval from —a to 
a. The circular sine (for which u-=1 and a==12/2) are simple 
illustrations of this case. 

In the case of the Jacobian function, in which 


Lt my +ayp=(l—yy(l—e*y) (0< k < 1), 
we have 
| 
a=. Q = | Se 
, Faw) d— ky) 
This last quantity is called the complete elliptic integral of 
the first kind in Legendre’s normal form, It is denoted by K (a) 
or, more briefly, by «A. Special tables have been compiled for 


it. We note that there is a simple expression for » in terms of 
K(1/V 2): 
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Thus, if we define the functions sin #f, sinh ¢, sl ft, and sn (f, k) 
as the functions inverse to the integrals referred to above, only 
the function sinh ¢ is defined on the entire real axis. The func~ 
tions sinf, sl ?f, and sn (ft, k) are defined only in the finite 
intervals 


w a + w . a a 
esis be. ga Koi tes. K, 


respectively. 

In the following chapter, we shall show how we can adhere to 
our basic plan of exposition in terms of the inverses of integrals 
and still extend these definitions not only to arbitrary real 
intervals but to the field of complex values of ft. However, when 
we do this, the product formulas and the periodicity of the 
functions sin ¢ and sl ¢t, of which we spoke in Chapter 1, will not 
be immediately obvious. Nonetheless, an additiontheorem, which 
we shall obtain in Chapter 4, will enable us to establish easily 
not only these but many other facts that go far beyond what we 
might discover if we confined ourselves to a geometric defini- 
tion of these sines. 


CHAPTER 3 


Integration in the Complex Plane 


12. As we have stated, it is possible to study completely the 
properties of the various sines and their similarities and 
differences if we treat them as functions of a complex variable 
to. 6--it, Where o and t are realnumbers and / is the imaginary 
unit. Of course, when we do this we cannot then interpret / as 
area or arc length, nor can we interpret the corresponding 
value of the sine in question as the ordinate of a point or the 
length of a chord. At the basis of the definition is the integral 


the upper limit w of which is some complex number ww. «# - dy, 
where « and wv are real numbers. 

Let us assume, just as in Chapter 2, that m and » are real 
numbers such that wm’ ---4a2+0. The integral (37) is a particular 
case of an integral of the form 


a 


[RG Warr ae raat are acts (38) 


where the integrand is a rational function of z and the square 
root of a polynomial of no higher than fourth degree. Such an 
integral is called, in general, an elliptic integral (provided the 
integrand does not reduce to a rational function of 2 and the 
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Square root of a polynomial of second or lower degree). In 
particular, for the integral (37) to be elliptic, it is necessary 
that m’—4n+40 and that 240. The name ‘‘elliptic’’ is due to the 
fact that the length of an arc of an ellipse is expressed in terms 
of an integral of the form (38). Inparticular, if »/a? + y?/o = 1 
(where (+a) is the equation of a particular ellipse, the length 
of the arc CA (see Fig. 13) is represented by the integral 


oe i Vi+(2yax ‘n= | i Baste ee 


Teh nae (a> — Av-) 


a 


1 i at — {a* —- b-) x” 
= nnn In 
V@ = x) (9 (a | 


FIG. 13. 


13. Let us pause to see how we should interpret the integral 
(37) in the case in which wis a complex number and what prop- 
erties this integral possesses. Here, we shall make certain 
assertions without proof, Proofs can be found in most texthooks 
on the theory of functions of a complex variable. 

Let us set up a Cartesian coordinate system in the plane. 
We shall regard (x.y) as the image of the complex number 
z==x + /y. In what follows, we shall apply the term ‘‘point’’ 
both to a point in the complex plane and to the complex number 
corresponding to it, 

Let us connect the coordinate origin z= 0 with any point 
w-==u iv by drawing a continuous curve L of finite length from 
the origin to the point w (see Fig. 14). (The fact that the curve 
starts at the coordinate origin is of no significance for the 
derivations that we are going to make.) Let us partition / into 
arcs by choosing points z,==0, 2. 2), ..., 2,.; %,=w numbered 
in the direction of motion along £ from 0 to w, and let us 
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form the corresponding approximating sum for the function 
1/V 1 mz? fanz: 
(39) 


pe er ey 


FIG. 14. 


Since V1 + mz’? +nz' is a double-valued function, let us agree 
that, for z= 2,0, we take the positive square root of 1. Sup- 
pose also that, if the value of the square root has been chosen 
at any point z,_,, that value at the point z, is chosen that will be 
the closer to the value at z,_,; Then, for an arbitrary partition 
of the curve L, the sum (389) will have a definite value. Let us 
now take over finer partitions of / so that the length of the 
longest arc approaches zero, It can be shown that the sum (39) 
will approach a definite limit. It is called the integral of 
1/V1 4 mz? nz‘ along L and is denoted by 


/ VI-+ m2? -+ mz oo) 


If the point # lies on the real x-axis (that is, if w= « and 
v=-0), if the polynomial | + mz*+-nz' does not vanish anywhere 
on the interval [9.4], and if L coincides with this interval, then 
all the values of z, that is, 2,0, 2,,..., z,, are real numbers, 
the values of Yl -+ mz" +-nzi are real and positive, and we have 
simply the usual definite integral 


1] 
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In analogy with the general case, it is customary to write the 
integral (40) in the form 


dz 
a —", 41 
| Vi+ mz? + nz =) 


without indicating the curve ZL over which the integration is 
taken (that is, for the points of which the approximating sums 
(39) are formed). However, such notation causes some in- 
definiteness. Specifically, for two distinct curves L, and L, that 
unite the points 0 and w, the corresponding integrals (40) may 
have distinct values. However, in the study of the theory of 
functions of a complex variable, it is shown that this can happen 
only when the curves L, and L, encircle one or more zeros of the 
polynomial | + mz*-+ nz‘ (at which the integrand becomes infinite) 
One Or more times. [See Fig. 15, where f, and p, represent two 
of the four zones of that polynomial and the curves L, and L, 
together encircle both these zeros twice.] Here, the difference 
between the two values of the integral depends on just how many 
zeros are encircled and how many times they are encircled. 
Therefore, the integral (41) has infinitely many distinct values 
for every value of w; that is, it represents a multiple-valued 
function F(w), any one of the values of which we shall denote by 
the letter ¢ (where ¢ is a complex number): 


t= / 2 = Fe). (42) 


FIG. 15. 


Let us suppose that w, does not coincide with any zero of the 
polynomial | + mw? nw’, that is, 
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1 + mus + nw} +0. 


Then, we can draw a circle X around w, that will also encircle 
no zero of that polynomial (see Fig. 16). Let us choose any 
curve L, connecting z=-0 with w, For points w inside the circle 
K, let us take curves formed by the union of L, and the line 
segment drawn from w,.tow. Under these supplementary condi- 
tions, the integral (42) defines a single-valued function of -w 
inside the circle, which coincides with some single-valued 
branch of the function; F(w) inside the circle kK. Its value at the 
point w is equal to the integral along L, plus the integral along 
the segment ww. To switch to any other branch, we need only 
replace the curve l) with another curve Ly to which another 
value of the integral at the point w, corresponds. It is easy to 
see that the corresponding branches differ by a constant term 
inside AK, One can show that each of these branches possesses 
a derivative at an arbitrary point w inside the circle K; spe- 
cifically, 


(43) 


FIG. 16. 


In other words, the theorem that the derivative with respect 
to the upper limit of integration of an integral is equal to the 
value of the integrand at the corresponding point remains valid 
for the complex integral (42). 

Furthermore, one can show that none of the values of our 
integral that are assumed at any point w, inthe plane can 
coincide with any of the values of the integral assumed at 
another point w, (where w,#w,). In other words, to a given 
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value f¢ of the integral there corresponds only one point was 
the upper limit of integration,* 

All this justifies our regarding w informula (42) as a single- 
valued function of the complex variable ¢, That is, we may 
regard the upper limit of integration as a function of the value 
of the integral. Let us denote this function by 


: w= S(t). (44) 


As before, we shall call it a generalized sine. 
It follows from (43) that 


dw =: ) 1-+-mw?-+nw' dt. (45) 


With regard to the multiple-valuecdness of the integral (42), 
which results from the fact that to a single value of w there 
corresponds infinitely many distinct values of ¢, this property 
leads to the following assertion for the inverse function (44). 
A single value of w is obtained for infinitely many distinct 
values of t. In connection with this, we shall see in Chapter 6 
that the function s(f) is a periodic function of ¢. 

14. Let us compare the values of the integral (42) corre- 
sponding to two points wand — wsymmetric about the coordinate 
origin. Let us connect these two points withthe origin by curves 
L and L’ that are symmetric about the origin (see Fig. 17), and 
let us integrate over these curves from the origin to the points 
wand --w, respectively. If we partition the curve L into arcs 
by means ot the points z,==0, 2, z)...., 2, wand evaluate the 
corresponding approximating sum 


v I (46) 


a eae end (Z, — 24. \)s 


and if we partition the curve ZL’ into symmetric arcs with the 


aid ofthesynmimetric points ~- 2,=:4. Big! Te se Rae o_o it 
and again evaluate the approximating sum, 


ee ee (46’) 


*A detailed proof of this fact can be found in FE. Coursut,-\ Course in Mathematical 
Analysis, vol. Il, New York, Dover, 1959. However, this proof assumes familiarity 
with the general theory of analytic functions and the fundamentals of the theory of el- 
lipric functions. 
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the values of the second sum willdifferin sign from the preced- 
ing one. But the first of these has the limit 


| daz 


b l-p mz" +n! 


and the second has the limit 


It is clear from what has been saidthat these values differ from 
each other only in sign. In other words, it 


then 


| ciel waite aes =5 t e 
x Pl peamz pact 
(We are still assuming that the integration is over two curves 
symmetric with respect to each other about the point > -» 9.) 
From this it follows that, if the value of the upper limit w 
corresponds to the value of the integral /, then the upper limit 
corresponds to the value -—/. In other words, if 


& fi) cw 
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then 
s(—)=—w, 
that is, 
s(— f= — s(P). (47) 


This means that the function s(¢), representing the inverse of 
the integral (42), is an odd function. 


FIG. 17. 


15. In the preceding section, we compared the values of the 
integral along the curves / and L’, where L’ was obtained from 
L by a rotation through the angle x about the coordinate origin. 
Now, let us compare the values of the integrals along the curves 
L and Ll”, where L” is obtained from L by rotating it through 
an angle «a/2 about the coordinate origin (see Fig. 18). We note 
that this rotation corresponds to multiplication of each point of 
the curve Ll (that is, the complex number represented by each 
point) by the number i. 
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Obviously, if the curve / connects the origin with the point 
w, then L” connects the origin with the point iw. 

Let us partition the curve JL into arcs separated at the 
points z) = 0, Z 2 -..+ Zaye 2%, = w. For the function 
1/V Tf mz? + nz’, the corresponding approximating sum is 


n 


1 
(Ze : 48 
= Vitm2+nz! Ve aie) 


The points obtained from the above-listed points by means 
of a rotation through x/2 about the ele eke origin lie on 
the curve lL”, These are the points iz,)=0, iz), iz. ..., iz, 4 
iz,=iw. The corresponding approximating are for the same 
function 1/V 1+ m= + mz is of the form 


n 


»  {LSp = £25.04). (49) 
a \, l— mt aes 


If the partition of the curve / is varied in such a way that the 
length of the longest arc z,_,z, approaches 0, the first approx- 
imating sum will approach the limit 


i ae tt 
; V i+ mz? + nzi ' 


which represents one of the values of the multiple-valued 
function 


/ Vi+ mz? + nzi : 


Comparison of (49) and (48) shows that (49) can be regarded 
as the approximating sum constructed for the same curve but 
for a different function, namely, the function 


i 


a fh 
Vi— mp ++ nz, 


Therefore, the limit (49) is equal to the integral 


Berea dz 
Vi—mz? nz — mz? -t-nzi’ 
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which represents one of the values of the multiple-valued func- 
tion 


w 
/ f dz 
3 Vi—- m2? -++ m2s i 


On the other hand, this same limit must be equal to the 
integral 
[ dz 


me V 14 mz? + nz 


which is one of the values of 


G 


since the approximating sum (49) was constructed for the curve 
L” and the function 


l 
Vit mz: + nz * 
From what has been said, it follows that for every value of the 
integral 


tw 
dz 
‘ V1+ mz? + nzt 


there corresponds the equal value of the integral 
aw 
if dz 
7 V1— mz? + nz? ° 
In other words, if 


iw w 
; V 1+ mz? + mz Q V 1 — mz? + nz 


then ¢=//f* under the hypothesis that the integration in the first 
integral is over the curve L” obtained from the curve L, over 
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which the second integral is taken, by rotating L through an 


angle x/2 about the coordinate origin. 
Let us denote s*(f) the function inverse tothe integral 


w 
f dz 
u V l— mz? + nz . 


Obviously, s*(4 is also a generalized sine. It follows from the 
equation 


w 
i} 2 UE = SS 
VY l—m2?+ nzé 
that w*=s*(#*). Also, iw—s(t)and t=if. Therefore, 


s(f)=s (it*)= lw = is*(t*). 
Finally, if we drop the asterisk in the notation #, we have 
a (if) = is’ (£). (50) 


This relationship connects the values of two generalized 
sines, the first representing the inverse of the integral 


/ dz 
J VT m2? + az | 


Q 


and the second representing the integral 


w 
f dz 
a Vi— mz? + nz! : 


16. Using these facts, in what follows we shall consider 
each of the special cases of the generalized sine as a function 
of a complex variable. Specifically, the circular sine w= sint 
is the function inverse to the integral 
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the hyperbolic sine w = sinh? is the function inverse to the 
integral 


w 
=f dz, 
ge Viz! 


the lemniscate sine w=s!f is the function inverse to the inte- 
gral 


w 
6 — 


the Jacobian elliptic function (the sine amplitude) w=sn(t,) is 
the function inverse to the integral 


w 
dz 
= ———[SSSSSSSS=—== 0 . 
lp Se 


If ¢ assumes real values, we can assert that the functions 
sin f, sinh ¢, and sl ¢ defined here coincide with the functions of 
Chapter 1 only in those intervals in which these last functions 
are inverses of the corresponding integrals. For sin ¢, this 
means that —n/2 <¢< n/2 and, for sl/, that —w/2<t< o/2. 
Only for sinh ¢ are we without any such restriction; in this 
case, the integral has an inverse for the entire real axis (cf. 
Chapter 2). 

From section 14, it follows that each of these functions is 
odd: 


sin(— ft) = — sinf, sinh (—t)=—sinht, 
sl(— f)==— sl f, sn(— f)—= -— sn. 


(51) 


To use the derivations of section 15, we need to consider the 
function s(t) together with the function s*(f), The change from 
s(f) to s*(¢#) corresponds to change of sign of the coefficient of 
z* in the polynomial | + mz?+ nz! in the integrand. To see this 
clearly, we present the following table of type sf) and s*(é) 
functions corresponding to the four particular cases that we are 
considering. 
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Integral Inverse | Integral Inverse 
{ 
ww w 
[ eee eee S$ (f) | ieee SAE s* (t) 
V1 — mz? +124 Vi—mz? + nz‘ 
w w 
| aoe 2. sin f | wince sinht 
| Viz Vi-+2? 
0 0 
w Ww 
| as sh é / a sin’ 
a } + « 0 l—2z 
Ww Ww 
: dz ‘ dz 
— — slr | slf 
| VP Yss2! Vi—2 
0 
ww va) 
F d F A . 
| ae sn (i. k) | fet ns sn (Z. k) 
; Vd -= 2°) (1 — #2?) . aa + 27) (1+ kz?) 


I 
| 


If we apply to each of these particular cases the general formula 
(50), we find 


sin(it)=isinhé, (52) 
sinh (it)=7 sin é, (53) 
sttit) = ist t, (54) 

sn (it, k)==isn*(é, R). (55) 


It is easy to see that formulas (51) follow from these last 
formulas, For example, from formulas (52) and (53) applied 
successively, we obtain 


sin (—t)=sin{i(it] = isinh(it) =1- isint=—sint. 


17. The results of sections 14 and 15 can be obtained ina 
very simple manner if we assume that the familiar rule of 
change of variable under the integral sign is also applicable to 
integrals of functions of a complex variable. In courses on the 
theory of functions of a complex variable, it is shown that this 
rule is indeed applicable under quite broad conditions—conditions 
that are satisfied by the examples we have been considering. 
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In the integral 


w 
° 


ral 
+'| 


| dz 
| Vt mz?- 


let us make the change z=—&, We obtain 
w ~w 


=| dz = , qe 
J Vipme fast Pian cae 


Turning to the inverse functions, we conclude that 


w—=—s(t), —wes(—t, ie, s(--f = --s(h. 


If we make a different substitution, namely, z = -~-/, that is, 
$= iz, we obtain 


Cel 
xT > . 
2 2 i dt 
t= — j ee 


When we turn to the inverse functions, we conclude that 


w= s(O, IW = SUD. 16.4 SG) = 18): 


Let us use once more the rule for change of variable, this time 
to express the function sn‘(f, k) representing the inverse of the 
integral 


0 
| TSU ere 


in terms of the Jacobian function sn¢ with appropriate modulus, 
To dothis, we make the following change of variable z= ¢/V1—@. 
After some simplifications, we obtain ‘ 


w bold we 
ie 2 


dt 
Vd-@T—d— ye 1 


| see ae ees 
1 V2) (1+ 2") 
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Let us set k’=Y1—k?>0 (where &’ is called the modulus 
complementary to k), Then, the equation 


_— —— 
f VIF 2 OF B24) 


which indicates that w= sn‘(f, k), can be rewritten 


fz f bien aan 
6 (0) (LG) 


from which we get 


w 
a ean (fs RB” 
Vi-w ) 
or 
VT sn? (tk) 


Therefore, equation (55) can be replaced with the equation 


Shura ae (55’) 


where &’ = V1 -- &. 


CHAPTER 4 


Euler’s Method of Deriving the 
Addition Theorems 


18. At the basis of all trigonometry lie the two addition 
theorems 


sin(a-+ fp) = sinacosf-+cosasinf, (56) 

cos (a-+- 6) = cosacosf — sinasinB. 

With the aid of the relationship between the sine and cosine, 

we can rewrite these addition theorems in such a way that the 

value of sin(a+ ) is expressed only in terms of sina and sinp 

and the value of cos(a-++- 8) is expressed only in terms of cosa and 
cos fp: 


sin(a + 8) =sina Vi — sin? p+ Vi- sin? a sin fp, 


Beant AP ope Se (56’) 
cos (a+ f) = cos acos p — Vi — cos’ a Vi — cos? . 


We do not write the + sign in front of the radical, but it should 
be understood that the sign is taken accordingto the value of the 
circular function expressed by that root. For example, 


V1 — sin’ = = cos Mas VS. 


6 2 
f 3x 3 V2 
—sjn? — =- eee 
V | — sin’ =- == cos == 5 ete. 


Each of formulas (56’) expresses algebraically the value of the 
function of the sum of two numbers in terms in the values of 
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that same function of the two individual numbers. Formulas of 
this type are called algebraic addition theorems and we say that 
a function for which such a theorem is applicable possesses or 
obeys an algebraic addition theorem. Thus, the sine and cosine 
obey algebraic addition theorems. Other examples are the linear 
function y:=:at and the exponential function y=:e', These are 
examples since 


a(a-+ B) =: aa -+ af, et 8 — 9%. of. 


19. Our next problem will be to derive an addition theorem 
for the function w= s(¢) inverse to the integral 


Ww 
(= (/ ——“— 
3 V14 mz? 4 nz! ; 


In particular, we shall obtain addition theorems for sin f, 
sinh ¢, sl vt, and sn(¢t, &) In this section, we shall observe the 
general behavior of the solution of this problem. 

Suppose that a and f are arbitrary complex numbers and let 
us define y=a-+p. We set 


u=s(a), v=s(p) w=s(a+f). 


This means that 


a= f dy p= | dy 

v Vip my ay! 2 Vi- my? anys” 
2! (57) 
rr 


y= f- 
| Vimy? fay’ 


0 


If a and # vary in such a way that their sum y remains constant: 


a +fp=y=const, (58) 
then 
d(a +p) =0, 
that is, 


iu 


° dy dy 
A Vit my? + ny! eS VI my? + ny" 
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or 


du. a dv = 
V1+ mu? + nu‘ Vist mv? + nv 


(59) 


This is a differential equation connecting the values uw = s(a) and 
v=s(f) subject to the condition (58). We shall show below that 
equation (59) has an algebraic solution and we shall find it. More 
precisely, we shall find an algebraic function F(u, v) that is 
equal to «w at v=0 (and equal to v atu=—O), the differential of 
which vanishes by virtue of condition (59) (or (58)). In other 
words, we Shall show that (58) implies that 


dF(u, v)==0 and F(u, v) =C =const. (60) 


To find the constant C, we set p—0. Then v=—s(f) vanishes and 
F(u, v) coincides with u=s(a). But, if B—0, it follows on the 
basis of (58) thata —y and hence u=s(a)=s(y), so that 


CF (a, Y) [geo = 4 |pey = SC) 
Thus, 
F(u, 0) = F{s(a), s(B)] = 5 (y) == 5s (a-+ f). 


This relation is independent of the value of the sum a+-fp=y. 
Consequently, we shall obtain an algebraic addition theorem for 
the function s(f): 


s(a—+ B) = F[s (a), 5 (B)]. 


Thus, the problem of deriving an algebraic addition theorem for 
the function s(f) reduces to finding an integral of equation (59) 
possessing the properties listed. 

20. To find an algebraic solution of equation (59), we shall 
use a method proposed by Euler in his Institutiones Calculi 
Integralis, Vol. I (1768).* We shall begin with the following 
fourth-degree algebraic equation relating the variables u and v 
and containing three arbitrary parameters A, B, and C: 


u? +- uy? 4- Aurv? + 2Buv -- C? = 0. (61) 


*Here, we shall use a simplified form of Euler’s method that will be sufficient 
for our purposes. 
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If we differentiate this equation, we obtain 
(u + Bu + Aunv?) du 4-(v + Bu + Au*v) dv = 0. (62) 
But, if we rewrite (61) in the form 
(Av? + 1) a? + 2Buu + (v? —- C2) = 0, 


multiply both sides by Av?+ 1, and complete the square, we 
obtain 


{(Av? + 1) u + Bol]? — [(C? -- v2) (Av? + 1) Bv2} = 0, 


from which we get 
a+ Bu + Auv? = Y C? + (B?+ AC? — 1) v2? — Avt. (61') 


Again, we do not write the + sign but we understand which- 
ever of the two values of the square root coincides with the 
expression on the left. (An analogous remark will apply also to 
the functions that we shall consider below when we encounter 
square roots.) 

Equation (61) is symmetric in « and v. Therefore, if we 
reverse the roles of # and v, we obtain from (617) 


vt Bu-+ Au?v = VCP (BP AC? — I) a? — Aut (61”) 
If we substitute (61’) and (61”) into (62), we obtain 


VOR PLAC ta dota + 
+ VC?+ (B?-+- AC? — 1) ua? — Autcv =0 


or 


ec a a een nS 
V C?-+ (B+ AC? — Tu? — Au! 
dv (59”) 

= 0. 

ag V CP -+ (B? + AC? — 1) ot — Avt 


Since the differential equation (59’) is satisfied for all values of 
«a and w that are related by the algebraic relation (61), it follows 
that (61) is an algebraic solution of (59%). 
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Let us now choose A, B, and C so that equation (59’) will 

coincide with equation (59) of the preceding section. To do this, 
we only need to set 


Bt AC2— 1—=mC2, A=—AaC. 


Let us express A and B in terms of m,n,and C. Then, 
equation (59’) (after we multiply all terms by C) takes the form 


du du 


Viimeta Vera (59) 


and its solution (61) is 


ue -+- U2 —nC2n2v? + 2 V1 + mC?+4+ 1Cluv —C2=0, (63) 


where C is an arbitrary constant. 

Since we wish to represent (63) in the form (60) (that is, in 
the form F(u, v)=C), let us solve equation (63) for C. We 
obtain successively 


{u2+ v? — (nav? + 1)C?2]2? = 4(1 + mC? + aC) u29?, 
(1 — nu2v?)? C4 — 2 [(1 + nu?v?) (u? + v2) + 
+ 2mu2v?| C2 + (a? — v?)? =0, 
(1 + nu?v’) (u? 4- uv?) + 2mu?v? 
irae = ~ (lesan 7 
x V {0 - nu?v*) (a? £0?) + 2mu?v? |? — (@? — 07)? (1 — nu? vy? = 
(1 — nu?v*)? 
_ uw? (1 + mv? + nv‘) + 0? (1 + mu? + nu‘) is 
a (1 -— nu?v?)? 


fi Quv V0 + mu? + nu’) (1 + mv? + nv) 


( — nuvy? 


so that, finally, 


uV 1+ mv? avi 4 ou Vd me? $n —C 
ee ea ee ae (64) 


The function 


Pa, 0) = AVE ma Ee + o VTE mi ak 


1 — nu?y? 


is algebraic. At v==0, it becomes u and, at w--0, it becomes vu. 
We leave the reader to verify that its differential is of the form 
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d du 
dF (a, v)== O(a, %)( u ), 


V i+ mu? + nus az V i+ mv? + anv! 


where @M(u, v) is some nonzero algebraic function. Therefore, 
the differential equation (59) implies 


dF(u, v)=0, or F(a, v)=const. 


Thus, the existence of an algebraic solution of equation (59) 
possessing the properties indicated in section 19 is established. 
From this, as was shown in section 19, we get an algebraic 
addition theorem for s(t) in the form 


s(a+ p)= F[s(a@). s(B)}. 
that is, 


s (a) V 1 ms? (B)--ns* (B) +s (B) VU ms? (a)-F asta) 
1 — ns? (a) s? (B) = 
s? (a) — s? (B) 


“OVE OO OVitwOTEe (6) 


s(a +B) = 


21. From the general addition theorem (65), which is valid 
for the function s(f), we obtain as special cases the addition 
theorems for sinf, sinh ¢, sl ft, and sn(¢, &) Corresponding to 
the circular sine are the values m -—— |!andn=—0; corresponding 
to the hyperbolic sine are the values m= 1 anda :=0; corre- 
sponding to the lemniscate sine are the values m=-0 and 
n-==— 1; corresponding to the Jacobian function are the values 


Accordingly, from formulas (65), we obtain the following 
addition theorems: 
For the circular sine, 


sin(a-+ B)= sine 1 — sin? f + sing V1 — sina; (65’) 


for the hyperbolic sine, 


sinh (a -+ 8) = sinha V1 -4- sinh? 8 + sinhB V1 -+ sinh? a: (65”) 


for the lemniscate sine, 


. slaVI— sh p+slRV¥T—sha 
he 1-Eslfast?p — 


a sPa—sl? fp . 
sla WI —si*B—sIPVI—sl*a ’ 


(65’”) 
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for the function sn(ft, &), 


sn(a + §)= 
__ snaVi-+ msn? p+ asn'p+sn fh V1 + msn?a--asnta = 


_ l—nsn?asn?p 


sn? a — sn? B 


— See ER Set PO ; IV 
snaV1I+msn’?B +asn'f —snBV1+msn?a+-nsn'a eer) 


For brevity, we write sn(a-+) instead of sn(a-+-f§. &) and sna 
instead of sn(a. k), etc. 

It is important to note that a and ff in these formulas are 
arbitvary complex numbers except that the a and # in formulas 
(65”) and (65'%) are such that the denominators of the fractions 
are nonzero, 


CHAPTER 5 


Further Study of Complex Values 
of the Argument 


22. For real values of ¢, we have the formulas 
cost== Y1—sin¢t, cosht = V1 -+ sinh’?, (66) 


where corresponding to the value f= 0 is the value of the 
radical equal to 1. Thus, formulas (65’) and (65’) take the 
forms 


sin(a +B) =sinacosfi-+ sinfcosa, (67) 
sinh (a -+- B) = sinha cosh8-+ sinh B cosha. (68) 


From these and from formula (66), we derive the addition 
theorems for the circular and hyperbolic cosines. Specifically, 


cos(a-+ 6) = V1 —sin’(a+ p) = 
= Vsin’ asin? Bp — 2 sinasinBcos «cos p + cos? a@ cos? B = 
= + (sinasinB — cosa cos f). 


To choose the proper sign, we need in this case only to set 
p=0 We see that we need to take the minus sign on the 
right: 


cos (a -+- B) = cosacos f — sinasinf. (69) 
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Analogously, for cosh(a 4-6), we obtain, to begin with, the 
formula 


cosh (a -+- 8) = Vi + sinh? (a + 8) = 
= V(cosh? a — sinh? a) (cosh? 8 — sinh? 8) + (sinha coshB + cosha sinh B)? = 
= Veosh? a cosh? 8 + 2cosha cosh sinha sinh 8 + sinh? a sinh? B= 

= +(cosha coshf + sinha sinh #). 
For p~—0, we obtain cosh a on the left and + cosh on the right. 
Therefore, on the right, we need to choose the plus sign: 


cosh (a + 8) = cosha cosh f ++ sinha sinh f. (70) 


Returning to formulas (66), let us replace ¢ with /f. We 
obtain 


cos (it)= V1 — sin? (id= V1 -+ sinh? t = cosht, 
(71) 
ete 
cosh (it) = V1-+ sinh? (it)= V1 — sin?t = cosht, 
(In the right-hand members, we choose the plus sign since, for 
t=-0, the values of both circular and hyperbolic cosines are 1.) 


If follows from formula (66) that both cosines are even 
functions: 


cos(-- t)== WY 1l—sin?t==cost, cosh (—t)= Vl -+ sinh?t=cosht. (72) 


From the addition theorem (67), we also conclude that the 
formula cos t=sin(x/2—f) is valid for all complex values of 
t, To see this, let us set a = x/2 and B= —¢ in (67). We obtain 


sin (5 — ¢) = 1 - cos(—#)+-sin(—#)- 0 = cost. (73) 
Here, if we replace ¢ with x/2 — ¢, we obtain 


COs (S$ —e)=sin t. " (74) 


Thus, formulas (5) of section 1 are valid (for all values of ¢ 
both real and imaginary) if we define sin ¢ as the function 
inverse to the integral 
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w 
— [ we 
§ Vl—z 


and define cos ¢ as the square root of (1 —sin?¢) (with the con- 
vention that cos 0 = 1). Since, under such a definition sin ¢ is an 
odd function [see formulas (51)| and cos ¢ is an even function 
[see formulas (72)], formulas (4) of section 1 also remain valid. 
But, as was shown in section 1, from formulas (4) and (5), we 
get all the ‘‘reduction’’ formulas of section 1, in particular, the 
periodicity of the circular functions. Therefore, in what follows, 
we shall use the formulas of section 1 for arbitrary complex 
values of f: 


im nr : 
cos ($+ 4) = ~sin t, sin (F +t) =coss; 


cos (m—~- t) = — cos f, sin(a — 4) =sinf; 
cos (x-+ f)==— cos, sin(m- ¢) = —sinf; 
cos (2 -} f) == cos f, sin(2n--¢) =sinf. 


To summarize, we may say that the new definitions that we 
have made (beginning in Chapter 2) for sin ¢ and cos ¢ have led, 
for real values of ¢t, to the functions that were defined on the 
entire real axis in section 1 on the basis of geometric con- 
siderations. We note that, for the hyperbolic functions sinh ft 
and cosh ?t, we do not need to go through the analogous proof. 
This is true because the difference between the cases of the 
circular and hyperbolic functions consists, with regard to the 
present matter, in the fact that sin ¢t, treated as a function of 
the real variable inverse to the integral 


2 
p=) es ae 
5 Vi-7 


was defined directly only on the finite interval —n/2<t< 12/2 
and we had no information, to begin with, regarding its behavior 
outside that interval. The hyperbolic sine, on the other hand, 
treated as the inverse of the integral 


¥y 
i= f gies, 
| Vite 
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was defined from the very beginning on the entire real axis 
(see Chapter 2). 

If we set a=o and §=/t in formulas (67)-(70), where o and 
t are real numbers, and use (52), (53), and (71), we obtain 


sin(a-+ i7)=sino cosh7 +i sinhr-coso, 


cos ( ir)=coso coshr—isino- sinhr, 
oe (75) 


sinh (o-++ i7)=sinho cos7-+-isin7-cosho, 


cosh(o + i7)=cosho cosr-+isinho: sinz, 


These formulas express the circular and hyperbolic func- 
tions for an arbitrary complex value of the argument t=o + it 
in terms of these functions of the real variables o and 1. 

23. In section 8, we derived expressions for cosh fand sinh ¢ 
for real ¢ in terms of the exponential function: 


t -¢ t “¢ 
e e - ewe 
ant, aa sinh t= 


cosht = 


From this it follows that 


e’ = cosht-+-sinh/. (76) 


Since we defined cosh ¢ and sinh ¢ above the arbitrary com- 
plex t—o-++dét, formula (76) enables us to define the exponential 
function for an arbitrary complex value of ¢. Using this de- 
finition, we obtain, with the aid of formulas (75), 


e7*!7 — cosh(a + in) + fsh(o+ in) = 
=cosr(cosha+ sinho) + isinr(cosha-+ sinho). 


Now cosho + sinho=e”° inaccordance with formula (76). Therefore, 


ef = eF +! — eF (cost + /sint). (77) 


If we replace ¢ with /f’/in (76), we obtain 
eit” cosh(it’)-+ sinh (it) = cost + isint’, 
or, Omitting the primes, 


e —cost+isint. (78) 
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This formula is due to Euler. Here, if we replace ¢ with —f, 
we obtain 


e—## — cost — tsint. (79) 
From (78) and (79), we derive two other formulas of Euler: 


it it C2. ad 
cost SA sint = (80) 


2t 
Let us verify that the addition theorem 


et+B — e2. @B, (81) 


where a and 6 are arbitrary complex numbers, holds for the 
exponential function e! defined by formula (76). We have 


e**8 — coshla + 8) +- sinh(a + B) = 
=(cosha - cosh8 + sinha - sinh8) + sinha - coshf 4+- sinh - cosha) = 
= (cosha ++ sinha) (cosh B + sinh B) = e?- e?. 


[Here, we used formulas (70), (68), and once again (76).] In 
particular, 1=e°=e#-4—e%-*, From this it follows that the 
value of e“ is always nonzero no matter what the value of the 
complex number «a. 

24. Let us turn to the lemniscate cosine cl ¢. In section 5, 
this function was defined for real values by the formula 


gees ($—~<). (82) 


Let us take (82) as the definition of cl ¢ for arbitrary complex ¢. 
If we set o=o/2 and p=—t in (65”) and note that sl o/2=1 
{formula (18’)] and si(— ¢)=—slf[formula (51)], we obtain 


_ ) exh ree [—sl?t 
ct=st(5—f)= ee pe ESS ar a (83) 


The value of the radical is chosen in such a way that cl O= 1, 
Equation (83) can be written in the symmetric form 


sPf+cPf+sPt-cPt=1. (84) 
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From this, we get 


l1—cl?f¢ / 
We note that, if we replace ¢ with it, we obtain from (83) with 
the aid of (54) 


wie, Len ey BSS 
en V ern Tae ae’ (eo) 


Furthermore, if we replace an arbitrary complex number 
fin (83) with —¢, we get, by virtue of (51), 


foes Ee sy l—sP?r 


Thus, cl ¢ is an even function. With the aid of (83), the expres-~ 
sions (65”) take the form 


sl(a+ p) = +25 B ( fare ahs a(1fsl?a) _ 


Pee RUE Las Gee ee (87) 


To establish the addition theorem for the lemniscate cosine 
as simply as possible, let us first replace ¢ in (82) with »/2~—1. 
We obtain 


cl (3 a t) = ee (82’) 
From (82), it follows that 


cl (a-+ 8) = si (5 —a] = p| 5 l($ eas B) be a]. 


Let us replace a by w/2—a and 6 by ~—fin the first expres- 
sions (87) and let us replace a by o/2— 8 and B by —a in the 
second expression. If we then use formulas (82) and (82’), we 
obtain 


lactB(1-+sl? Bp) —slasl B(1 4+ cl? 
clu ¢ By = SER Se 
cP? B—sl?a 


™ Clacip(-+sl2a)+slasifp(l elf) * (Se) 
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In contrast with the expressions (87), these last fractions 
are not represented symmetrically with respect to a and f. 
However, the formulas are of course valid when we reverse the 
roles of ao and 6. 

25. We have already noted that all the ‘‘reduction formulas’? 
of the circular functions are consequences of formulas (4) 
expressing the evenness of the function cos ¢ and the oddness 
of the function sin ¢ and of formulas (5) expressing the equality 
of the values of cos ¢, and sin ¢, if t, +t, 2/2. The lemniscate 
formulas possess analogous properties: cl? is an even and 
sl f is an odd function [cf. formulas (86) and (51)]. Further- 
more, the values of cl ¢, and sl ¢, are equal if ¢,+4,=— w/2 (cf. 
formulas (82) and (82’)]. Therefore, the corresponding ‘‘reduc- 
tion formulas’’ obtained by the natural substitution of sl for 
sin, cl for cos, and o for 1 are valid: 


cl (3+ 4) = -- se. si($+t)=as 
cl(a — ft) == —clf, s!(w—-¢f)==sl ft; 
cl(a + f)ae—clf, SI (w-+ f) = -- slf; 

cl(2w-+ f) = cl é, sl (20 + f) = slit). 


It follows from this, in particular, that the lemniscate func- 
tions are periodic functions with period 2a. Now we may assert 
that, for real values of t, the lemniscate functions defined by 
considering the inverse of the integral 


y 


t= f il 
leea 


on the entire real axis coincide with the lemniscate functions 
introduced in section 5 by way of a geometrical approach, 

The identity (84) shows that, for real values of ¢, the values 
of the lemniscate functions 


x=clf and y=slt 


can be regarded as the coordinates of a point A lying on the 
following fourth-order curve (see Fig. 19): 


xP yP 4 x2y? = 1, (84’) 


Together with the point A (x, y), let us consider the point 
A’(é, y), where § =1/x. In accordance with (85) and (54), we have 
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l 1 1 


Obviously, — and y satisfy the equation obtained from (84’) by 
replacing x with 1/£ but leaving y as it was: 


(84) 


FIG. 19. 


We have thus obtained a new fourth-order curve (one con- 
sisting of two unbounded branches as shown in Fig. 19). The 
coordinates of an arbitrary point A on it are expressed in terms 
of the values of the lemniscate functions of purely imaginary 
argument it: 


SES. Ss 4 sl (é0). 


The variable ¢ has an easy geometric interpretation. Spe- 
cifically, from equations (84’) and (84”), we have 


[-y?@ ity? 


XS Sey = ———— 
Yi-y: VI—yi 


(For definiteness, all the points are taken with positive abscissa.) 
Therefore, for the areas of the curvilinear trapezoids OCAK 
and OCA‘K, we obtain 


y y 
: - 1a 
area OCAK= | xdy=> | —=—— a7, 
/ J i 
y y 
a ‘° 2 
area OCA'K = | oan == | we = 
U 0 oe 
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from which we get 
pe fee NES ct * [area OCAK + area OCA K]. 


Such a geometric interpretation of ¢ in this method of repre- 
senting the lemniscate functions is valid both for the case of a 
real parameter ¢ and a purely imaginary parameter it. 

The addition theorems enable us to express the values of the 
lemniscate functions for an arbitrary complex variable ¢ =o + it 
in terms of the same functions of the real variable o and 1. 
Specifically, when we set ao and £ = it in the first of expres- 
sions (87) and use the formulas 


sl (it) =Ts! 7, (54) 
(85) 


: 1 
cl (it) Saar 
we obtain 


slo (1 — sl? t) + ési rclo (i -+ sf? a) 
1—sl*o sl? t = 7 


si (o 4+ it) = 


In this equation, let us replace ! —sf°t withcl’t(1--sl*t)in ac- 
cordance with formula (84). After dividing through by cl t, we 
obtain 


slocl t (1 4+-sl? 1) -- fsl tcla(1 + sf? 9) (89) 
Th sla Slee ° 


si(o + ft) = 


Analogously, the second of expressions (87) yields, for a =o and 
B= it, 


so -b sl? t 


/ 
SO eH) are ch peat a) el Clo size) * oe) 


In a similar manner, we derive from formulas (88) the relation 


cloclt(l+sl?t)—ésloslt(l+cl?o) _ 
L—cPosl? x a 


l—sPocl’t 
= Samay I ee) 


cl(o + it) = 
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26. Let us now turn to the formulas for the Jacobian elliptic 
functions. The basic function sn(f, &) is used to define two other 
Jacobian elliptic functions: 


cn(t, k)== V l—-sn?(t, k). dn(t. k)=z WI—R? suet, &). (91) 


Here, we take values of the radicals that assume the value 
1 at ¢--0, that is, at which sn(0, &£)-= 9, and we extend the defini- 
tion of these functions to the complex plane in such a way that 
continuity and differentiability are ensured. Thus, 


cn (0, k) == dn(0, &) == 1. 
The relations (91) can also be represented in the forms 
sn?(¢, k)--cn2(4, Ry=— 1, kesn2(t, ky) + an2(t, RY 1. (91’) 


It follows from formulas (91) that cn(¢, k) and dn(f, &) are 
even functions 


cn(— f, k)==cn(f, k). = dn(—?t, k)= dnt, ). (92) 


The notations for the Jacobian functions sn ¢ and cn ¢ under- 
score their similarity with the functions sin fand cos f, respec- 
tively. If the modulus = 0, these functions reduce to sin ¢ and 
cos f. The function dn ¢ then reduces to the constant 1. 

To see this, note that for k =0, the integral 


aw 
dz 


f= eee 
| V (= 2?) (l= #2?) 
becomes 
bg 
pss | ae 


Therefore, the inverse function w=sn(t, k) becomes w—sinft, 
Formulas (91) yield ‘ 


en(f, 0)=cost, dn(f, 0)=1. 


With the aid of formula (91), we can represent formula (55’) 
in the form 
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sn (t, 2’) 


sn (éf, 1 aT RTS 


(93) 


where k’ = V1 —k? is the complementary modulus (with respect 
to £). Formulas (91) now enable us to express cn(ét, k) and 
dn(ift, k) in terms of the corresponding functions of the comple- 
mentary modulus: 


7 e«nt(iF. by tk l 
a, Da VIR ME Be V+ Se ara! OM) 


dn (it, A) = Wl— Fk snte, & =Vi+k ea = 


_ Vis (t, ’) + (1 — kh”) sn? (4. B) dn (£, &’) 
cn (t, k’) Yeh dey,” (95) 


If in formulas (93)-(95), we replace ¢ with éf (and hence zt 
with —¢), we obtain 


sn (it, k’) 
en (if, R’)’ 

1 __ dn (it, &’) 
coun ey? 99C ar a 


sn(f, k)=—l 


cn (f, k) = 


Noting that the complementary modulus k’ vanishes for k=], 
we obtain, with the aid of formulas (52) and (71), 


a2 sn(it,0) ss, sin(ét) __sinht | 
sat =~ bade 0) ‘cos it) ~ cosh? ° 
; _ Bi 2 et) ee. 
cn(f, 1) en (if. 0) cos (if) cosh (f) ° 

dn (it, 0) | ! 


date cn (it, 0) = Cos (if) ~~ cosht * 


Consequently, for k=1, the Jacobian functions reduce to 
ratios of hyperbolic functions. The corresponding graphs (for 
real values of f) are shown in Figure 20. Thus, the Jacobian 
functions snit, &), cn(t, &), and dn(t, k) for continuous variation 
in the modulus # from 0 to 1 form, as it were, a bridge from 
the circular functions sin ¢ and cos ¢ and also the function 
identically equal to the constant 1 to the hyperbolic functions: 
(sinh s)/cosh 1, 1/cosh ¢, and 1/cosh ¢, respectively. 
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en(t,1)=dn(t,1) = ——— 
ee 


FIG. 20. 


27. With the aid of the functions cn(¢, 2) and dn(¢#, &), the 
addition theorem (65!) for sn(f, &) takes the form 


su(a+-p, &) == 
7 sn (a, | (a, ek) cn (f. k) dn (B, &) + sn (B, &) cn (a. &) dn (a, ky 
~ } — k? sn? (a, k) sn? (fp. &) 


a sn? (a, &) — sn? (B, 2) (96) 
~~ sn (a, k) cn (PB. &) dn (B, &) —sn (B, &) cn (a, &) dn (a, &) | 


To derive addition theorems for cn(f, k)and dn(t. k)from (96), 
let us use the following identities, the validity of which can be 
checked by use of formulas (91’): 


(snacnfdnfp-—+ snfpcenadna)?+ 
+ (cnacnf —snasnf dnadnp)? = 
= (dnadnf — k?snasnBcnacnp) + 
+ k(snacnBdnfp-+ snfhcna dna)? =(1 — &* sn’? asn?’ fp). 


We obtain 


cn(a +-f, AE) = VI— svr(a+f, k)= 

sa V0 — &? sn? asn’ py? — (snacnp dnp + snfcna dna)? = 
OT Re sn? a sn? B Se 

= etacnp sagen fl dnednp., (97) 
_ 1— k? sn? asn?B 


dnfa £ Bf oat Ca k) -= 


_ Fa- - Rk? sn? a sn? Bp)? — &? (snacnBdn fp + snf cna dna)? 
= | — &? sn? a sn? B 4 


__ dnadnf —&?’ snasnBenacnp — (98) 
_ ] — k* sn? asn?p 
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These addition theorems enable us to express the values of 
the Jacobian elliptic functions for an arbitrary complex number 
é—o-+it in terms of the same functions of the real variables 
o and t+ Here, we also need formulas (93), (94), and (95). The 
first of equations (96) yields 


sn(o-+ /t, k) = 

__ sn (a, &) cn (ét, Rk) dn (it, &) 4+ 8n (it, &) cn (9, &) dn (G, Rk) 

= | — &? sn? (co, &) sn? (it, ) 

_ sno, &) dn (t, 2’) + isn (t, &’) en (1, 2’) cna, &) dn (a, A) (99) 
= cn’ (t, Rk’) + Rk? sn’ (co, k) sn? (t, &’) 


where k’ =: ¥1—&?. The second of expressions (96) enables us 
to represent the same result in a different form: 
sn(o + it, k) = 


sn? (a, 2) cn? (1, &’) 4+- sn? (t, 2’) (997) 
sn (o, &) dn (1. k’) —éisn(t, &’) cn (t, &’) cn (a, &) dn (a, &) 


In the same way, by using formula (97), we see that 


cn({o + it, k) = 
cn, k) cn (t, k’) —isn (a, &’) sn (t. &’) dn (o, &) dn (1, Ri) (100) 
cn? (t, k’) + 2 sn? (co. &) sn? (1, &‘) 


The reader can use formula (91) to verify that the same 
result can be represented in the form 


cn(o-+ it, k) == 


_ cn? (o, k) + &” sn? (o, &) sn? (1, &’) (100’) 
~" en (a. k) cn (t. k’)-+ isn(a, &) sn (t, k’) dn(o. &) dn (t, 27) © 


Analogously, for the function dn ?, we obtain 


dn(o + it, k) = 

__ dn(o, &) cn (t, &’) dn (t, k’) —- ik’ sn (9, &) cn (a, &) sn CE, Ry 

oF cn? (1. Rk’) + &? sn? c k) sn? (t, 2’) 

= ke cn? (a. k) ER” cn? (1, &’) (101) 
~“dn(o, Aen (a, &’) dn (t, &’) + ck? sn(o, &) cn (o, &)sn(t, &’)* 


28. Let us set 


pees! = K(k) =K, (102) 
Vl — y4) 0 — Fy?) 


where the integration is along the segment of the real axis 
between the points 0 and 1, so that K(&)is a positive real num- 
ber. It follows from (102) that 


sn{K(k), kJ =|, (103) 
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Therefore, 


cn[K(k), A] = V1 — sn’ [K (Rk), kA] = 9 


and 


dn [K (k), R} = V1 —k?sn?[K(k), AJ = VI — B=” (104) 


For simplicity of writing, let us drop the symbol for the 


modulus & If we then seta—=K andfp = —f in formulas (96)-(98), 
we obtain 

sn K cnédnt—sn¢cnk dnk cnfdnt cn ft 

sn (K — t) = 1—&? sn? K sn? ¢ = To#R suit ~ dnt’ (105) 
, cnAcnt4+-snK sntdnK dnt 
SL Se) eaalnaea Nou orca eT) aaa 
Rk’ sntdnt , sine | 

japan ar. -e2) 

_ dnKdnt+ksnkKsntcnKenf  k’ dnt _ k’ 
a 1— k? sn? K sn? t ~ Jk sn't dnt’ (107) 


If we now replace f by —f, we obtain 


sn(K +f) =<2o = sn(K Set), | 

, sunt ¥ 
no ag era (108) 
dn (K -- j= s == dn(K — f) 


In the theory of Jacobian elliptic functions, formulas (105)- 
(108) play a role analogous to the reduction formulas in trig- 
onometry. 

If we replace ¢ by f+ XK in formulas (108), we find further 
‘‘reduction formulas’?: 


snif-+ 2K)=—snf, | 
cn(f + 2K)-=—cnf, (109) 
dn(f-+-2K)=—dnt. | 


c 


From the last formula, it follows that the function dn fis 
periodic with period 2k, Finally, in the first two of formulas 
(109), let us replace ¢ by t-+ 2k, This gives us 


sn(f+4K)==snf, cn(f-+4K)==cnt. (110) 
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Consequently, sn ¢ and cn ¢ are also periodic functions—in 
their cases, with period 4k, 

29. The formulas of the preceding section enable us to 
characterize completely the behavior of the Jacobian functions 
of a real variable without any other calculation. Let us begin 
with the function y=sn(t, kh) (for0<&< 1). From the definition 
of this function, specifically, as the function inverse to the 
integral 


dy 
: ernest 


0 


it follows that it increases monotonically from 0 tolast 
increases from 0 to 


1 
= en 
KW= | 7S 


(cf. section 11). The first of formulas (108) shows that the 
graph y=sn(t, k) is Symmetric about the straight linet—K 
and hence has the shape of a half-wave in the interval 0<t< 2K 
(see Fig, 21a). Furthermore, if we replace ¢ by —-¢ in the first 
of formulas (109), we get 


sn (2K — t) == — sn(—f)=sné. 
By comparing sn(2K +-¢) with sn(2K — 15), we obtain 
sn (2K + t)== — sn(2K — 0). 


This means that the graph of y=sn(t, k) is symmetric about 
the point (2K, 0). Consequently, in the interval 0-<¢t< 4k, it has 
the same shape as does a sine wave inthe interval 0<t<2n 
(see Fig. 21a). Since we know that the function sn(é, 2) is 
periodic with period 4K, we can now proceed to construct its 
graph. 

To determine the nature of the graph of the function v= cnif, k), 
we can begin with the relation cn(f, k)== Yl —- sn(t, &), from 
which it follows that this function decreases monotonically 
from 1 to 0 in the interval 0 <1< A. Furthermore, the graph 
is constructed just as above, by use of formulas (108)-(110). 
We see that this graph looks much like that of the cosine (see 
Fig. 25b). 
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dn(tA) 


FIG. 21, 


In an analogous manner, we obtain the graph of the function 
y=dn(t, k). The formula 


dn(t, A) = VI —ksn°s, &) 


shows first that this function decreases from 1 to Y1—k#? =k’ 
in the interval 0<t<Kk. The graph is then extended to the 
entire real axis by using the appropriate relations (108) and 
(109). It is shown in Figure 21c., 

It follows from the above analysis that all real zeros of 
sn(t, &) [that is, the real roots of the equation sn(¢, k)-=-0] are 
exhausted by points of the form t=2mK(k), where mis an 
arbitrary integer, 

Analogously, the real zeros of cn(f. k) are of the form 


fon (2m — 1)K(R). 


The function dn(f, &) has no real zeros at all Gf 0<k< 1). 

30. Let us show that the generalized sine s(¢) can in all 
cases be represented in terms of the functions sin ¢ and sinh ¢ 
or in terms of the Jacobian functions, Specifically, let us show 
that, for 7-=0 (where n is the coefficient of z‘1 in the expression 


FURTHER STUDY OF COMPLEX VALUES OF THE ARGUMENT 65 
1+ mz?-+-nz*), it can be expressed in terms of sin?/or sinht 
and that, for n # 0, it can be expressed in terms of the Jacobian 
functions. 


1+22z2, where A>0. Then, if we represent the relation 
w— s(t) in the form 


tw 
r= | dz 
od V1 + A2z? ‘ 


and make the change of variable €-=22, we obtain 


>| 


Aw 
1 dé 
f= =—, 
| Vi + 


From this we have either }w =sin(2.t) or Aw= sinh (A/), Thus, in 
the case in question, 
s()=+sin(at), or s(t) = ~ sinh(A, t). (111) 
Suppose now that n+0 and, as always, that m’?’—4n2+0, We 
have two cases: 
CASE I. m?— 4n>0. Here, the roots x, and x, of the equa- 
tion 


x? -+ mx +-n == 0 
are real and distinct. We have 
x? mx mse (KX — Xj) (KX), 
or, if we replace x with 1/z-, 
Lo mz?*taztis=() - x ,27)(1 52°). 

At this point we need to consider three subcases: (a) x, and x. 
both positive, (b) x, and x, of opposite sign, (c) x, and +, both 
negative, 

Subcase (a): We can assume that x, and x, have the values 


7x2 and 2, where 2? > > 0, Then, 


1+ mz? nz4 == (1 — 2.227) (1 — p?z?). 
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Consequently, the equation w—=s(t) takes the form 


2 ae 
t= | -. = ; 
| VU — 222770 -— pez?) 


If we make the change of variable 6= Az, we obtain 


Aw 


feat f aM 1 syhere0 Sb ve lL, 
Ae =e} a 


Thus, in this subcase, 
w= s(t) == sn(ht, R). (112) 


Subcase (b): We can assume that x, = 4? andx, = — jt’, where 
4 and | are both positive. We obtainthe following decomposition 
into factors: 


Let. mz? + nat=(l- xy222Q1 — xyz?) a (las) (1 + we’). 


If, in the integral 


ess | ot dz 
‘ Vd — 2222) Cl 4 22?) 


we make the change of variable C= \1 -- 22z?, we obtain 


l Dcpee? 
| ,* 1 
a | Wo = = 
pa pe V2) = RY) 
where 
(0) ko eiose Paso ae | 
VAY fe 
If we set 
,] 
"ee — 


e FOS) re) 
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we obtain 
Vi -22my? 


K- V A? 4 pet de 


Vd d=") 
or 

VT hie = sn(k— VE = 
__ n(Va? +p? t, k) 


in subcase (b), so that 


mV? Ppt, ky ke sn (WRU tk) 


l 
woe § (f) == — 1 -- = om eS ; 
©) rs dn PA? + pit ki i dn(VA24p?e, &) 


Noting that, in the present case, 


we can rewrite our last result in the form 


! sna? yet. kd 


cml cma cceerea Tt cera he at) 
For 4 =p = 1, the integral 
@ 
Dida ee ot ee 
v VL aez*) (1 + p27) 
coincides with 
f de 
Vie 
Therefore, s(f) = s1(t), and we obtain 
a 1 
sil (¥ fee aol 
sl ft . — Lise ie (113’) 
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For cl tf, we obtain 


dn? > 
ata y/ it mt (V 28s 
pst = on (V2 


= cn (V2. 75) (113”) 


Thus, the lemniscate functions can be expressed in terms 
of the Jacobian functions with modulus 1/V2. 

Subcase (c): When x, and x, are both negative, they can be 
represented in the forms —i? and — py’, and we can assume that 
A >p>O0. Therefore, 


1+ mz*+-nz4=(1 — xz?) (1— 927) = (1 4 A*z*) (1 + 22). 


In the integral 


_ | : dz 
dh VF 22?) (2?) 


we make the change of variable € = 1z/V'1 +222? (cf. calculation 
in section 17, where 4-= 1 andjt=—&). We then obtain 


dt Vip? 
(= ideas 2_______where0 < & = ———- <_ |. 
A i) Vd =e (1 20?) A > 
From this, we get 
Aw ‘ ‘ 
Viper = Sn (Af, k), 
or 
iene l sn (At, 2) __ | sn (At, k) (114) 


A V1 —sn? (Af, 2) A sn (At, R) 
It remains for us to consider 


CASE Il. m? —4n <0 (obviously, here n > 0). Since the roots 
x, and x, of the equation 


7+ mx t-n==0 
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are imaginary, the roots of the equation 
£4 + m§?-++n2=0, 


defined by 


2 $2 __ 
G=G=—*, &=& 


are also imaginary. Furthermore, their real parts are nonzero 
(because otherwise the squares of the roots would be real num- 
bers). Noting that ¢,, &, & and & must also be pairwise conju- 


gate and changing the numbering if necessary, we have 
SiS Ke +h B= —-&=—a—h @>0,p> 0). 


In the decomposition 


ES + mE n= ~ ENE — BYE — ENE — 80) 
let us replace £ with 1/z. We then obtain the following decom- 


position into factors: 


I-+ mz? + nzt== 
=|1 + (a + if) z] (1 +4a—ép) 2] {1 —(a +48) z] [1 —-(a—@p) zJ= 
=[I 4 2az + (a? +p) 27] (1 - 2az + (oe? + B) 27}. 


Obviously, a? + p?= Wn (from which we getO<a< Wn) Con- 
sequently, our integral is of the form 


dz 
z?)(1 -Qae + Van 2?) 


te 
pap 
6 V (1 4-2az + V1 
Let us show that it reduces to one of the cases considered 


above when we make the substitutior 


When we make this substitution and perform some manipulations 


in the integral, we obtain 
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70 
4 
14+Vaw 
4 
VYaw 
i ‘ | ac = 
= | "a 4 ta \ ee (i= vo: 
: (veal Va qa)ee] | Vn +u}p Va—a)e| 
4 
leh nw 
a 
lib aw 
odie { de — 
p . V+ 262) (1 2) 
fa ee 
pea SE en Ee 26 
; Va —u Vn +a 
Let us set 
| 
f ae mS i, ane 
“ Vl ER) T+ pe?) 
Then, 
earn 
ayy we 
Bo te ee A 8 
THREAT EWE 


pPr+A= | 
r 
Thus, this brings us to Case I (c). Therefore, in accordance 


| sn (B+ ALA) 


with formula (114), 
4 
Aa Vnw 
4 en (Bi + AL RY’ 
l--foanw 
k= Vw 
z “ 


where 


sn (pt -+ A, k) —Acn (Bt + A, &) 


From this we get, finally, 
sn (ft + A, k) + Acn(pt+ A.) * 
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Thus, we have shown in this section thatthe generalized sine 
with 2+0 can be expressed in terms of the Jacobian functions 
with appropriate modulus & and, for 2-0, interms of the circu- 
lar or hyperbolic sine. All this will give us the right in what 
follows to confine ourselves to a study of only three sines 
(together with the corresponding cosines): the circular, trig- 
onometric, and Jacobian (the sine amplitude). However, together 
with these, we shall also consider the lemniscate functions since 
it is for this special case that the properties of the Jacobian 
functions manifest themselves in the simplest manner, 


CHAPTER 6 


Zeros and Poles. 
Simple and Double Periodicity. 
The Concept of an Elliptic Function 


31. If we compare formulas (75), obtained above for circu- 
lar and hyperbolic functions with formulas (89), (89’), and (90) 
for the lemniscate functions or with formulas (99)-(101) for the 
Jacobian functions, we see a sharp difference between formulas 
(75) and the other six. Specifically, the circular and hyperbolic 
functions are represented as polynomial functions of sin o, 
cos o, sinh o, and cosh ‘co, whereas the lemniscate functions were 
represented as fractional expressions in terms of sl o, clo, 
sl t, and cl t and the Jacobian functions were expressed as 
fractional expressions in terms of the Jacobian functions of the 
real variables o and t 

Therefore, the functions sin /, cos t, sinh ¢, and cosh ¢ are 
defined and have finite values for every complex ¢. On the other 
hand, the functions slf, clf. sn(¢, &£), cn(t. &) and dni¢, &) are not 
defined for those values of ¢ at which the denominator of the 
fraction in question vanishes, 

Let us look into this phenomenon more closely, beginning 
with sl¢, The denominator of the fraction in formula (89) is 
equal to! — sl*o-sl’t It vanishes if and only if 


shoc $lLt— + 1, (115) 


But slo and sit are real numbers that do not exceed 1 in abso- 
lute value. (We recall from section 5 that slo is the signed 
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length of the chord of the lemniscate.) Therefore, equation (115) 
is equivalent to the two equations 


shoes 2) and sl t= 1, (116) 


It follows from the same geometric representation of slo that 
the condition s!g= 1 is equivalent to the condition 6= w/2 + 2mo, 
where m is an arbitrary integer (2 being the period of slo). 
Since slo is an odd function, the condition 0 =—1 is equivalent 
to the condition o—= — wo/2+ 2mw. Thus, the denominator of the 
fraction (89) vanishes at all points with coordinates 


and only at such points. 
The corresponding complex numbers f=o-+ it are exhausted 
by the formula 


t= (im tS 4424 NS. 


But if m ranges over all integers 0, +1, +2, +3, ..., the quan- 
tity 4m + 1 assumes all odd values +1, £3, £5,.... 

Therefore, the values of ¢ can finally be represented in the 
form 


i 
t= (2p — 1) 54+ (2y— 1) (117) 


where p==0, £1. £2,... andg=0, +1, £2,.... 

Let us look at one of these points, If o== (2p — 1)w/2, then 
slo == (=a as can easily be derived from observations made 
above. Therefore, sl?o=1, and formula (83) yields clo=0, 
Analogously, for t—(2qg—1)/2, we have slt==(—1)7"" and 
clt==0,. If we substitute these values into the right-hand side 
of formula (89), not only the denominator but also the numerator 
in the formula will obviously vanishand we cannot draw any con- 
clusion as to the behavior of sift as t—>(2p— 1) w/2 4-(2¢ — L)wi/2 
immediately. Formula (89’) saves the situation. Here, the 
denominator of the fraction is equal to 0, just as before, but the 
numerator si?o+sl?t=2 and thus is nonzero. From this it 
follows that sl¢ approaches >> as ¢->(2p —1)0/2+(2g—1) wi/2. 
(Here, we use the fact that the lemniscate functions of a real 
variable are continuous.) We can now extend the definition of 
silt, to the entire complex plane by setting 


ZEROS AND POLES 75 


sl [ep — DF +9 — DF] =. 
PaO. Sol ee es ge Os ET. Oe es (118) 


If we refer to a point in the complex plane at which a func- 
tion becomes infinite (that is, at which the function approaches 
oo. as the argument approaches that point) as a pole of the 
function, we may now say that formula (118) provides all the 
infinitely many poles of the lemniscate sine. 

Analogously, it would be possible to determine all the poles 
of the lemniscate cosine from formulas (90). However, it is 
simpler to start with formulas (82): 


Ww ‘ 
cit==si(F—¢}. 


It follows from this formula that clf approaches co as ¢ ap- 
proaches a point a if and only if sl(w/2—¢) approaches ©» as 
»/2—ft approaches w/2—a. From this, we get, just as above, 

( 


5a = 2p —YVS +029 — DSi 


that is, 
a==2(—ptdS+12—9 + D—N-z- 


If we replace — p+ 1in this equation withp’ and if we replace 
—g+1 with q’, we see that all the poles of cl¢ are represented 
in the form 


t= p’w-+ (2q’ — D>. (119) 
where p’ and q’ are any integers. 

32. By an analogous procedure, we find the poles of the 
Jacobian functions. Let us begin withsn(¢, k) (where0<& < 1). 
The denominator of formula (99) is of the form 

cn?(t, k’)-+ &? sn? (a, &) sn? (1, k’). (120) 

Since the values of the Jacobian functions are real when o 

and +t are real, vanishing of this sum is equivalent to the two 


conditions 


cn(t, k')=0, sn(o, &)sn(t, k’) = 0. 
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But if cn(t, k’)=0, then sn(t, k’) = +1 [cf. (91)]. Therefore, we 
see that (120) vanishes if and only if the two equations 


cn(t, k')==0, sn(o, &)=0. (121) 


hold simultaneously. From this it follows in accordance with 
section 29 that 


o= 2mK(k) and t=(2n — 1) K(k’); 
that is, 
f= 2mK (k) + i(2n — 1)K (k’). 


Denoting K(k) by K and K(k’) by Kk’, we rewrite these values in 
the form 


t= 2mK +i(2n ~1)K’, (122) 


where m and a are any integers. 

However, if we substitute the values that we have found for 
o and t into formula (99), we see that not only the denominator 
but also the numerator vanishes [since sn(o, &)=cn(t, k’)=O}. 
Therefore, this formula does not enable us to draw any con- 
clusions immediately regarding the behavior of sn(¢. k) as 
t—>2mK + (2n —1)ik’. Formula (99’) comes to our aid. Since 
sn(t, k’)== £1 when cn(t, k’)==0 [by the first of formulas (91)], 
we have sn*(t. &’)=1, and formula (99’) shows that sn(t, &)—>0o 
as ¢ approaches any of the points (122). Consequently, the 
values found for ¢ represent the poles (all of them) of the 
Jacobian function sn (t, &). 

By an analogous procedure, we can see that the values (122) 
constitute all the poles of the functions cn(¢, 2) and dn(fé, k), 
Furthermore, the fact that the poles of these functions must 
coincide with the poles of sn(¢t, &) follows from identities (91). 

33. In what follows, we shall also need to know the zeros 
of the functions that we are studying, that is, all the points in 
the complex plane at which the functions vanish. Let us begin 
with the circular sine. From the first of formulas (75), we have 


sint = sin(o + ir) = sinocoshT + isinh7 coso. (75’) 
From this it follows that the equation 


sint=0 (123) 
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is equivalent to the system of equations 


sing - coshr.= 0, 
sinh; - coso = 0. (124) 


Let us look at the first of equations (124). Since cosht +0 
(recall the geometric representation of cosht in section 2 or 
formula (32): cosh t=e'+e-‘)/2), we conclude that sino~0, 
Therefore, o=—mn, where m=0, +1, +2,.... Consequently, 
coso =(—1)” #0. The second of equations (124) yields sinh 1 —0, 
so that t=0. (Recall the geometric representation of sinh t in 
section 2). Formula (31) leads to the same conclusion. Thus, 


t=o+it=mn, m=0, £1, £2,..., (125) 


This expression takes care of all the zeros of sin t. We see 
that none of them is imaginary. Extension of the definition of 
the circular sine to the complex plane does not yield any new 
zeros. 

We leave it as an exercise for the reader to use the second 
of formulas (75) to show that all the zeros of cos t are repre- 
sented by the expression 


=ofit=(2m—1)F, m=0, $1, £2..... (126) 
Finally, by the same procedure, we can find the zeros of the 
hyperbolic functions. However, it is simpler to achieve this end 
by using the relations 

sin (it) = isinht, cos (it) = cosht. 

It follows from the first of these that the equation sinh t—0 is 
equivalent to the equation sin(it)=0. Therefore, it—mn and 
hence 


f= (— m)uv—nni, 


where n is also an arbitrary integer. In just the same way, we 
can show that the zeros of cosh ¢ are given by the formula 


t=(2n—1)=, 


where nis an arbitrary integer. 
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34. Let us now seek the zeros of the lemniscate sine. Using 
formula (89’), we conclude that the equation si ft=s!(o + it)=0 
implies that sl?o+s!*1=0. Thus, we obtain the two equations 


sla==0, 
sl t= 0. (127) 


From this we get o-- mo. t=aw (where m and na are arbi- 
trary integers). Consequently, 


t=mo-+nol, m=0, +1, £2,...; n=0, £1, £2,..., (128) 


However, for slo=0 and sit=0, the denominator of the 
fraction (89’) vanishes. Therefore, formula (89’) does not tell 
us the values of sl¢ at the points that we have found. Let us turn 
to formula (89) for si(a-+/t), Here, the numerator of the fraction 
is equal to 0 when equations (127) are satisfied but the denomi- 
nator is equal to 1. Therefore, sl(mm-{nwi)=-0 and, hence, 
formula (128) does indeed give the zeros of si¢ [all of them 
since formula (89’) implies that the conditions found are nec- 
essary for vanishing of s!¢]. 

To find the zeros of cif, the simplest procedure is to use 
formula (82): 


a) 
Les (5 — 1). 
Cc S [= 


From this formula, it follows that equation clt-=0 is equivalent 
to the equation sl(o#/2—1)=0. Therefore, o/2—t—=mo-nol. 
From this we get 


tm (2(— m+ 1) — I] 5 +(—ajol= (2m! — 1) 24 n'ol, (129) 


where m’ and n’ are arbitrary integers. Formula (129) provides 
all the zeros of clé. 

35. Formula (99’) is convenient for finding the zeros of 
sn(f, 2). From this formula, we see that a necessary condition 
for sn(o + it, k) to vanish is 


sn?(g, &)cn?(t, &’) + sn2(t, k’)=0. 


This condition is equivalent to the two equations 


sn(o, k)cn(t, k’)=0, 
sn(t, k’)=—0. 
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But the second of these, together with (91), implies that 
cn(t, k’) = t1#0. Therefore, the necessary conditions for 
sn(o -+-it)to vanish are 


sn(o, k)=0, 130 
sn(t, k’)=0. ( ) 
Therefore (cf. section 29), 
O = 2mK (k), T== 2nKk (k’); 
that is, 
f= 2mK + 2ntk’ (131) 


where m and a are arbitrary integers. 

However, formula (99’) is not suitable for evaluating sn(f, &) 
at these points since, when conditions (130) are satisfied, not 
only the numerator but also the denominator of the fraction 
representing sn(f, &) vanishes. Therefore, let us turn toformula 
(99). Here, the numerator is equal to zero but the denominator 
is nonzero [we recall that en(t, #’)== +1]. Thus, (131) represents 
the zeros of sn(t, &) (all of them). 

With the aid of formulas (100’) and (100), the reader can 
show that 


f= (2m — 1) K + 2nik’ (132) 


provides all the zeros of cn(ft. &), and he can use (101) to show 
that 


f= (2m —1)K 420 — Wik’ (133) 


yields all the zeros of dn(f, &) (where m and nr assume all inte- 
gral values). 

36. With regard to the periodicity of the functions that we 
are studying, up to now we have considered only their real 
periods. However, in extending the definition of the functions to 
the complex plane, imaginary periods can also arise. An exam- 
ple of this kind is the exponential function. 

This function ef as a function of a real variable is not 
periodic. However, it does have the purely imaginary period 
2n. To see this, let ¢=~o+t, where o and t are real numbers. 
Then, from formula (77), 


ef == eF+!t — eT (cost + /sint). 
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Let us replace ¢ with f+ 2m. Since ¢-+ 2ni—o-+ (t+ 2n)i, this is 
equivalent to replacing t with t-+ 2n, We obtain 


ef t2mt — gt t+(t +20) — 9% (cos (t + Qn) + ésin(t + 2n)] = 


= e (cost + ésint)= ef 


for arbitrary ¢. From this it follows that 2ni is a period of the 
function ef, 

Let us recall that a number A is called a period of a func- 
tion f(t) if the equation 


ft+A)=F 


is satisfied for every ¢ for which the function is defined. 

If f(f) has at least one nonzero period, it is said to bea 
periodic function. It follows from the definition that, if Aisa 
period of f(t), then any integral multiple of A is also a period 
of f(t). 

It was shown in section 22 that 2n is a period of the function 
sin f, treated as a function of a complex variable. Consequently, 
any integral multiple of 2x is also a period of that function. Let 
us show that the converse is true: Any period of sin ¢ is an 
integral multiple of 2n. Suppose that A is a period of sin 4, 
From the identity 


sin(¢ + A)=sinf 


we obtain, for r= 0, 
sin A =sin0 =O. 

Consequently, A is zero of the circular sine. According to 
formula (125), A must be equal to some multiple of 1; that is, 
A-=mn for some integer m. Let us denote by & the remainder 
obtained when we divide m by 2. Then, m=2p+k, where p is 


an integer and k is either zero or l. 
Therefore, 


sin(t + A)=sin(¢+4+ en + 2pm) =sin(t+ kn). 
The fact that A is a period of sin ¢ implies that 
sin(f + ka) =sint 


for arbitrary ¢. If k=-0, this condition is satisfied. If ek —|1, it 
cannot be satisfied identically since sin(¢+)=—sint, Finally, 
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we conclude that k=0, that is, that A=2pn, which was what we 
wished to show. 

In an analogous manner, we can show that any period of cos ¢ 
is an integral multiple of 2x. In summation we may assert that 
the transition from a real to a complex variable does not bring 
about any new periods in the study of circular functions. 

The situation is different when we study hyperbolic functions. 
Formulas (52) and (71) provide the quickest path to what we 
wish to show: 


sin(if) = isinh¢, cos (it) = cosht, 
or 


sinh(t) =— rsin(it), cosh(t) = cos (it). 


From these formulas, we have, for an arbitrary complex num- 
ber A, 


sinh(t +. A)=—isin@t + 1A), cosh(t -+ A)=cos(it + iA). 


Therefore, A is a period of the hyperbolic functions if and only 


if i1 is a period of the circular functions, that is, if iA - 2pn, 
where p is an integer or, in other words, if A =—(— p)2ni = 2nip’, 
where p’==0, #1, £2, +3, .... Consequently, the hyperbolic 


functions sinh ¢ and cosh ¢ are periodic functions, and all their 
periods are given by the formula A==2p’ni, where p’is an 
arbitrary integer. 

We see that the transition to a complex variable disclosed 
the periodicity of the hyperbolic functions (these are not periodic 
if we consider them as functions of a real variable). 

37. Although the periods 2p. of the circular functions are 
all real numbers and the periods 2p’ni of the hyperbolic func- 
tions are purely imaginary (except for the period 0), both these 
classes of functions have in common the fact that any period of 
the function in question is an integral multiple of some single 
number. In the case of the circular functions, this number is 
21 (or -—2x); in the case of the hyperbolic functions, it is 2ni 
(or -—2zxi). If every period of a periodic function is an integral 
multiple of some single period, knownas the fundamental period, 
the function is said to be a singly periodic function. Thus, the 
circular and hyperbolic functions are singly periodic functions. 

Let us turn to the lemniscate functions, beginning with the 
lemniscate sine. Noting that slw = 0 and 


cla = sl (5 — ») v= 8 (—>) == — gs) —~=: —l, 
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we obtain the first of expressions (87) 
sl (f a a) sz — $l ft. (134) 


We see that o is not a period of sif. From formula (134), 
it follows that 


sl (¢-+ 20) = st [(f + o) 4 of] = —- sl (f4-0)= sit. (135) 


Therefore, 2o is a period of the lemniscate sine (as we noted 
in section 25), and hence every integral multiple of 2w is also 
a period. However, periods of the form 2 do not exhaust all 
the periods of the lemniscate sine! To see this, let us note that 
sl(w)=islw=0 and cl(wi)==1l/cdo = —1. If we set «=f and 
f =i», in the first of formulas (87), we obtain 


sl (¢-+ wi) = — slit. (136) 


From this it follows that «w/ is not a period of s!*. However, 
we conclude from this same formula (136) that 


slid + ad) =. sl{(¢ + ow) + oi) == — sltl + wip == sl bt; (137) 


that is, 2i, and hence any integral multiple of 2wi, is a period 
of sit. Thus, the lemniscate sine possesses infinitely many 
real and infinitely many purely imaginary periods. From this 
it follows that slit, being a periodic function, is not a singly 
periodic function. To see this, let us suppose the opposite, 
namely, that each of its periods is an integral multiple of some 
nonzero fundamental period u. Then, there must exist nonzero 
integers m and nv such that 


20 == ma, Qwi === Na ‘* 


But this would then imply that /== n/m, that is, that zis a real 
number, which is not the case. Thus, st¢ belongs to the special 
class of periodic functions that are not singly periodic. 

By using formulas (134) and (136), we conclude that 


sl (t+ w + fw) = sl f{(f + im) +] = — sl (f+ Wo) =slf. 


Therefore w + /o is a periodic of sit, as in any integral multiple 
of w+ lo, 

38. From the periods that we have found, we can construct 
new. ones by composition, For example, from the fact that 2 mw 
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and n(w-+ iw) are both periods of sit (where m and 7 are both 
integers), it follows that 


2mw + n(w + wl) (138) 


is also a period of sl ¢. 

Let us show now that every period A of the function s!¢ can 
be represented in the form (138) for some integral values of 
m and n. Of course, for those periods that we have already 
found, this is obvious: 


20=2-1-04+0-(09+4 wi), (M1, n= 0); 
2of = 2-—1-m42-(0+4+ 0), (m-==—1, n= 2); 
@-+ wi = 2-0-0 4+ 1-(w+ wi), On == 0, m==1), 


However, our assertion is of a general character and it 
states that sl/ has no other periods than those that can be 
represented in the form (138). 

Suppose that we have the identity 

sl(f-+ A)==slt. (139) 
Setting ‘== 0 in this identity, we find 
slA=0, 


From this it follows that A is a zero of the function s!tf; that is 
[cf. (128)], 


A= mw + nwi, 


where m and 7 are integers, Let us rewrite A in the form 


A= (m — n)w-+-n(w-+ ol) = po-+a(o+- ol). 


Let & denote the remainder obtained on dividing p by 2; that is, 
p=2r-+k (where r is an integer and & is either 0 or 1). Then, 


A = 2ro + 2(@ 4+ wi)-+ ko, 


When we substitute this value into (139) and note that 2rw-+ 
+n(o@+ i) is a period of the form (138), we conclude that 


sl (¢-+ Ro) ==sIit 
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identically. This condition is satisfied by only one of the two 
possible values of k, namely, k-=0 [that k —1is impossible is 
shown by formula (134)]. 

Thus, any period A of the function s!¢is of the form 


A= 2ro+n(o+ a), 


as we wished to show, 

The above analysis has brought out the fact that an arbitrary 
period of sit can be represented as a linear combination, with 
integral coefficients, of the two fundamental periods 2m and 
w+wi. Periodic functions possessing such a property are 
called doubly periodic. Thus, sl(f) is a doubly periodic function 
with fundamental periods 2m and » + wi, The reader may verify 
that we could also have chosen w— wi and » + wi, for example, 
as our fundamental periods. 

By using the formula 


clt=sl (> - t), 
the reader can easily derive from what has been shown that cl? 
is also a doubly periodic function with the same fundamental 
periods. 
39. An analogous analysis can be made for the Jacobian 
functions. Let us consider, for example, sn(t, &). In section 28, 
it was shown that 


sn({f/+ 2kK)=—snfé, 


from which it follows that 


sn(f-+ 4K)== sf; 
that is, 4K is a real period ofsn/. Let us show that snr also has 


imaginary periods. From formula (131), which gives all the 
zeros of sni(f, &), we conclude that ¢-= 2ik’ is also a zero: 


sn(2iK’, k) = 0. (140) 


Furthermore, from formulas (94) and (95), we obtain, for t-= 2k’, 


OK", k’ 
on (QiK’, Ry = ener dn (2iK’, y= Oe. (141) 
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But from formulas (109), it follows that, for r= 0, 
cn(2K, &)==-—-cn0=— 1, dn(2K, kA) =dnO—l, 


or, if we replace & by &’ and hence 


1 1 
: dt . f dt 
kel yesmacrs =! yaar 
we have 
cn(2K’, k’)==— 1, dn(2K’, Rk’) = 1. 
Therefore, formulas (141) yield 
cn(2K’, y= — 1, dn (2iK’, k’)= I. (142) 


If we set a=—f and p= 2ix’ in the first of expressions (96) 
and use formulas (140) and (142), we obtain 


sn(f-+ 2iK’, k)==sn(t, k)cn(2iK’, k)dn(2ik’, ky)=sn(t, k). (143) 
From this it follows that 2/kK’ is a period of sn(¢, 2). From what 
has been said in the present section, it follows that the complex 
number 

A= 4Km +-2iK’n, (144) 
where m and »n are arbitrary integers, is a period of sn(f, 2). 
Let us show that numbers of the form (144) exhaust all the 
periods of sn(f, 4). Here, we can reason as in the case of the 
lemniscate sine (cf. section 37). Let A denote any period of 
sn(f, k). Then, 

sn(¢+ A)—snt (145) 
identically. If we set ¢=0 in (145), we obtain 

sn A—OQ; 


that is, 4 is a zero of the function sn¢t. Therefore [cf. (131)], 


A= 2mK + 2K'n, 
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where m and ” are integers. Let r denote the remainder ob- 
tained when we divide m by 2: 


m=2p-+r, 


where p is an integer and r is either 0 or 1. Then, 
A=4pK + 2nik’ + 2rkK. 


If we substitute this expression into (145) and note that 4pK -+ 2nik’ 
is a period, we see that 


sn(f¢-+ 27K) == snt. 


This relation cannot be satisfied identically if r—1 since 
sn(t-+2K).= -- snt. Therefore, r—0, so that 


A= 4pK - 2nik’, 


which was what we wished to show. - 

From what we have proven, it follows that sn(¢. &) is a doubly 
periodic function and we can take 4K and 2/k’ as the fundamental 
periods. Analogously, one can show that cn(ft, &) is a doubly 
periodic function with fundamental periods 4K and 2K -+ 2/k’ 
[another convenient choice of the fundamental periods of cn(f, &) 
are the numbers 2k — 2ikK’ and 2k + 2ik”] and that the function 
dn(t, k) is also a doubly periodic function with fundamental 
periods 2k and 4ik’. 

40. Let us represent on the complex plane all the periods 
of sin¢ and let us draw parallel lines through these points. These 
lines partition the plane into strips (see Fig. 22), known as 
period strips. Obviously, these strips fill the entire plane 
without any gaps and without any overlap. The inclination of 
these lines can be chosen arbitrarily. Let us make the simplest 
choice and draw them parallel] to the imaginary axis. If a point 
f is allowed to vary in such a way as to assume all values in 
one of the period strips, the point (+ 2px (where p is an integer) 
will assume all of the values in another period strip. Since 
sinf assumes the same values as ¢ and ¢+ 2pna, all the values 
assumed by the circular sine in one ofthese strips are repeated 
in any other strip. Therefore, for example, when we seek the 
roots of the equation 


Sint fA, (146) 
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where A=:.a-+-é8 (a and f real) is a given complex number, we 
need only find its roots belonging toany one of these strips. The 
other roots are obtained from those found by displacements of 
the form (¢’=t-+2pn (where p is an integer). They occupy in 
each of the strips the same positions as the original roots do 
in the original strip. 


ps 
-2n \0 


FIG. 22. 


Let us show that, for any A=:2-+/if, equation (146) has 
exactly two roots in each period strip. (These two points may 
merge into a single point, called a multiple root.) Here, we 
need to adopt a convention assigning the boundary line between 
two strips to one or the other of these strips. Thus, let us 
assign each such line to the strip lying to the right of it, so that 
the points of any strip satisfy the inequalities 


2up<a< 2n(p +1) (—o<ct< + oo). 


For example, the points ¢-—0 and t=-/ belong tothe strip 
0<0<2n, and the points 2n and 2n-+/ belong to the strip just 
to the right, namely, the strip 2n <0 < 4n. Representing sint 
according to the first of formulas (75), we rewrite equation 
(146) in the form 


sing - coshr-+icos- sinhr = a+ if. (146’) 

We wish to show that there exist only two roots f:-0-+ it 

satisfying the condition 0 <¢o-< 2x. Obviously, (146’) is equiva- 
lent to the system 


sing-coshr=a,  coso- sinhr=f. (146”) 


Let us set sino=x, sinht=y. Then, these equations take the 
forms x V1 + y?=aand VY! — x*y=£ or 
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x24 y2y2— G2, y2 — x2y? — p?, (147) 
We are seeking the real solutions of (146”), Therefore, 


O<x?<! and 0<y?, These conditions are satisfied only by the 
values 


ee Fe te ee il (148) 


where 


A =: (a? + p? — 1)?+ 4p? = (@? + p? + 1)? — 4a? (149) 


and YA denotes the positive value of the square root. 

It is obvious from (149) that the values that we have found 
are nonnegative. To show that x*< 1, consider the difference 
1 — x?: 


(VE +p 
: | 


It follows from (149) that this difference is also nonnegative. 
The equation 


y? = sinh? pa OTR I1+ VA 
2 


yields two values of t that differ only in sign. Corresponding to 
them are two values of opposite sign, sinh t and —sinh 1, of 
the hyperbolic sine but only a single value of the hyperbolic 
cosine. 
The equation 
x? = sin?o = es 
is satisfied by two values of sin o, However, as the first of 
equations (146”) shows, we should take only the one whose sign 


is the same as the sign of a (since cosht > 1 > 0). We denote 
this value of sin o (which is the only one possible) by x,: 


SING==Xy, [xX _| <}. 


Furthermore, corresponding to this value of sin 9 are two 
values of o in the interval 0<o < 2x, which we denote oa, and 95. 
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The corresponding values of coso, and coso, differ in sign: 
COS 0, = — COS Op. 

Since the product cos o sinh t must be equal to $, replace- 
ment of o, by o, in the second of equations (146”) must be ac- 
companied with the replacement of one value oft by its negative 
—t. Thus, we have seen that there exist exactly two pairs of 
values of o and 1 that satisfy (146”) or (146) under the additional 
condition that 0<o < 2n: 


Oo, t and o., —t. 
The corresponding points f are 
t;=0,+it, t;=0)— It. (150) 


Here, o, and o, are determined by the conditions 


snot / SEP INVS 
j= ae 


(for j==1, 2), the sign of sino, coincides with the sign of «, 
O<co, << 2a, and + +t is determined by the condition 


sinh7 = + yf otha 1+Va_ 


We shall not write out explicit formulas for o and rt although 
this would not be difficult. 

We note that the points ¢, and t, are symmetric about 1/2 if 
a > 0 and that they are symmetric about 3x/2 ifa< 0. 

An analogous analysis could be made for thecircular cosine. 
[The simplest procedure for doing so would be to use the fact 
that cos ¢ = sin (x/2—f], 

Let us observe briefly how the results are changed if, 
instead of the circular sine, we consider the hyperbolic sine. 
Here, all the periods are purely imaginary, except for 0: 


OG; ea, Ant 


Therefore, for our period strips, we take strips lying parallel 
to the real axis (see Fig. 23). The identity 


sin (it) =: isinht 
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implies that the equation 
ginht = A, t=o+ir, O0< r<2z, 
is equivalent to the equation 
sin (it) = sin(—- 1 + io) =sin[(2xn —1)+4 lo] = 1A, 


where 0< 0’ = 2n--1 < 2n, We have already examined this equa- 
tion. By using the result obtained for sint,, we see that the equa- 
tion sinh t= A, for an arbitrary complex number 4A, has two and 
Only two roots in each of the period strips of the hyperbolic 
sine. (These roots may merge into a single double root.) 


T 


41. Let us indicate on the complex plane all the periods of 
the lemniscate sine: 


2mw-+ n(w + iw), m and n integers. 


In our studies of the periods of the circular and hyperbolic 
sines in the preceding section, we partitioned the plane into 
period strips. Here, we are dealing with a doubly periodic 
function and, instead of period strips, we naturally have period 
parallelograms. One of the possible choices of a system of 
period parallelograms for si is shown in Figure 24. As a point 
f assumes all the values in one of these parallelograms, the 
point t42mw+n(o-+ in), where m and vn are any fixed integers, 
assumes all the values in another parallelogram. Since 


sl [f4-2mu-+2(w + iw)| = sé, 


the values of si¢ assumed in one of the period parallelograms 
are repeated in every other parallelogram. Therefore, for 
example, in seeking the roots of the equation 
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sif{=< A, (151) 


where A—a- ép is a given complex number, we need only find 
all the roots belonging to any one of these parallelograms, for 
example, the parallelogram with vertices 0, 20, 20-+(w-+ ww), 
and +i, All the other roots of equation (151) can be found 
from these by means of displacements of the form 


M=t+2mo-n (wtio) . 


In each parallelogram, these other roots occupy the same posi- 
tion as do the roots found in the original parallelogram. 


Let us show that, for any A~=a-+/f, there are two and only 
two roots in a single parallelogram (these two may merge into 
a single point, i.e., a double root). Here, we again need to adopt 
some convention regarding the boundary points that are common 
to two or more parallelograms. Let us assign to each period 
parallelogram those points lying on the lower base (other than 
the right end-point) and the points onthe left side (other than the 
upper end-point). 

In the parallelograms I, II, III, and ITV shown in Figure 24, 
for example, the point A belongs to I, the point 8 to II, the 
points C and G to III, and the points Dand E to IV. 

Returning to the theorem mentioned above, let us first prove 
the following proposition: 

The rectangle 0<0< 20, 0<1<w (see Fig. 25) always con- 
tains two and only two roots of equation (151). 
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The original formulation of our assertion will follow from 
this since the function s!¢ assumes, by virtue of its periodicity, 
the same values in the triangle OAS outside the parallelogram 
OBEC, as in the triangle CDE. 


Let us express sl¢=sl(o+/t) in accordance with formula 
(89) and replace (151) with the equivalent system 


slo cl t(1 + sF 1) Sar } 
l—slPosl? 1 eee | (152) 
sl tcla (1 + sl’ a) 


~~ T—sPosl?t =p.| 


Let us express the values of clo and clt, respectively, in 
terms of slo and s!t in accordance with formula (83) and let us 
set slo= x and s!lt=:y. Then, the system takes the form 


xpl-—yi ; 
re 
pests | (152’) 
yp l— x! 
“Tye oP 


We are seeking only real solutions of this system that satisfy 
the conditions |x| <1 and 0<y< 1 (since sit>>0 for0<t< w), 

Let us square each of equations (152’), add the result, divide 
both numerator and denominator of the sum obtained by 1 -- x*y’, 
and then divide the whole equation through by 1 — x?y3 This yields 


x24 y* = (a?-+-B7) (1 -- x72), ay? CL = ct) -— p22) — yt)=0. = (153) 
We note that this equation remains unchanged if we simultan- 


eously reverse the roles of x and y on theonehand, and @ and £ 
on the other. Eliminating «? from (153), we obtain 
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pey® — [ (a? +62)? — 1] y — 2 (2a?-+ p?) yt + 


+ [ (a? + p?)? — 1] y? + p? = 0. eH 
If we divide through by y4 and make the substitution 
yore 
we obtain the following quadratic equation for ua: 
pa? +-[ (a> + p*) — |] a — 4020. (155) 


Since the conditions of the problem admit only values of y 
that are real and do not exceed 1 in absolute value, the quantity 
u==1/y?—y? must be a nonnegative number, Therefore, we ob- 
tain from (155) the single admissible value of u: 


wes Vier Py Ai Fe le tp 3.0: (156) 


Now, to find y, we have the equation 
yt + uy? —1=—0. (157) 


Here again, Only one root, a positive one, is possible for y’. 
Obviously, this root is less than 2 in absolute value (since the 
sum of the two possible values of y’, which is equal to —ug, is 
nonpositive). Consequently, it will not exceed unity: 


O< y?= Uae sa #1, (158) 


Since y=-slt>0O (where 0<1t< 0), this gives us the only 
permissible value of y. 

Analogously, if we set 1/x?— x?=:v and use the symmetry 
(mentioned above) of the system (153), we obtain 


ei A oh leet (159) 


¥ 2a 


and hence the only permissible value of x?: 


‘a2 LA , 
O< x? = — rae (160) 
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From the value that we have found for y= slIt: 


/ 42 —4 — 
yy) Ws 


we find two values t, and t, of t in the interval O0<t<, These 
two values are symmetric about o/2: 


We note that the corresponding values of cl t, and clt, are nega~ 
tives of each other (cf. section 5). For definiteness, let us 
suppose that y-=slt,=slt,>0. (We suggest that the reader 
himself examine the case y==(.) Then, from the second equa- 
tion in the system (152), it follows that the sign of clo coincides 
with the sign of 6. Therefore, for any two possible values o, and 


o, such that 
Voqpa—v 
slo, = sl 0 = y/ Verein 


where o, > 0. a, > 0, 6, + 0, =o, we should choose that one (which 
we assume to be o,) such that the signs of clo,and Bp are the 
same. Analogously, of the two values o, and a, such that 


ee eee 
; Vvw+4—v 
sie, =sto,=—|/ ¥ eka, 


where @-<6,, <6, and o,-+ 0, == 3, we choose the one (which we 
assume to be o,) such that the signs of clo, and B are the same. 

From what was said in section 5, wecan see that o, -+ 0, = 20, 
(We begin with the fact that 0-< 6, <0, © <6; < 20, sing, = — Sindy, 
and clo, -=clo,.) 

To summarize, we have shown that the system (152) can be 
satisfied under the conditions 0<0< 20, 0<1t< by only two 
values t, and t, of t such that ‘ 


slt;=slt,=yv>0, O<clt,=—clt,(t, +1, =), 
and by two values oa, and o, of o such that 


slo, ==— slo, >0, clo,s=clo, (a, +0, = 2), 
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Here, the sign of clo, (and that of clo,) coincides with the sign 


Since the first equation of the system (152) implies that the 
Sign of slo-clt coincides with the sign of a, it follows that the 
Signs of sio and clt must be the same if a > 0 and opposite if 
a<. 0). This means that the value o, is compatible with only one 
of the two values 1, and t,, namely, with t, in the case a > 0 and 
with t, in the case a >0. Similarly, the value o, is compatible 
with t, or with 1,. 

What all this amounts to is that there are two and only two 
points ¢, and ¢,-in the rectangle 0<o< 20, 0<1< » that satisfy 
equation (151). In each of the possible cases, namely, 


t=O + lt, ty == 0) +4 it, = (20 — 0) + ((~—1,) (a > 9), 
tp2=0,- ity, fy =-0,+ i, =Qo0—4,)+l(o—,) (a < 0) 


the points ¢, and t,; are symmetric about the center w+ iw/2 of 
the rectangle. 

If both these points do not fall simultaneously in the parallelo- 
gram OBEC, then one of git let us say, ¢,, belongs to the tri- 
angle OAB and the other, ¢,, belongs to the parallelogram OBEC 
(Fig. 26). But, in this case, the point 4,4, + 2, falls in OBEC, 
and we again have two roots of equation (151) in that parallelo- 
gram. In this case, they aresymmetric about the point 20 + iw/2. 


A B D 2 
o a 
FIG. 26. 


42. An analogous analysis could be made for the Jacobian 
functions with arbitrary modulus. Here, we would need to show 
that, for arbitrary A, any equation of one of the following forms 


sn(f, R)==A, cn(f, k)== A, dn(t, A)—A 


has two and only two roots (which may merge into a single 
multiple root) in the corresponding period parallelogram. How- 
ever, the proof would require some tedious calculations, and we 
shall not go through it. 


96 THE REMARKABLE SINE FUNCTIONS 


We prefer to conclude our exposition with a brief survey of 
the definitions and general properties of elliptic functions. 
These functions constitute a further generalization of the func- 
tions that we have been studying. At the basis of this study is 
the concept of an entire function of a complex variable.* 

A function f(f) of a complex variable is called an entire 
function if it is defined and differentiable at every point of the 
complex plane. Examples of entire functions are all poly- 
nomials, the exponential function (e'), the circular functions sin¢ 
and cosf, etc. It can be shown that the class of entire functions 
coincides with the class of functions that can be represented by 
everywhere-convergent power series: 


f(t) == gt ayt+ ay? + ag... -agt™ +... 


A function g(t) of a complex variable is called a meromor- 
phic function (from the Greek words uépoo ‘‘fraction’’ and 
nopen ‘“‘form’’) if it can be represented as the quotient of two 
entire functions: 


_ git) 
si 


Examples of meromorphic functions are all rational func- 
tions (i.e., the quotient of two polynomials), the functions 
sec t=1/cos ¢, and tan ¢=(sin ¢)/cos ¢. It turns out that sl4¢, 
cif, sn(¢t, &), cn(f, &), and dn(f, 2) are also meromorphic functions, 
However, this assertion requires a special proof and does not 
follow immediately from the fact that these functions can be 
expressed by fractions of the forms (89), (90), (99), (100), and 
(101). The difficulty lies in the fact that the numerators and 
denominators of these fractions are not entire functions of the 
complex variable t=o-+-/t. The proof that the lemniscate and 
Jacobian functions can be represented as ratios of entire func- 
tions and ways of finding such entire functions are special 
questions that would take us beyond the scope of the present 
book. Of course, an entire function f(f) can be regarded as a 
special case of a meromorphic function since f (t)= f (f/1. 
However, such a function has a finite value at every value of 
t without exception, whereas a meromorphic function of the 
form o(t)= ¢(4)/h( has, in general, poles at which it becomes 
infinite. Obviously, every point ¢=—f, at which A(t))=0 but 
g(t)) # 0is a pole of the meromorphic function ¢ (4). 


“See, for example, the author’s Entire Functions, New York, American Elsevier 
Publishing Co., 1966. 
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A meromorphic doubly periodic function is called an elliptic 
junction The above implies that the lemniscate and Jacobian 
functions are elliptic. 

One of the basic theorems in the theory of elliptic functions 
asserts that every nonconstant elliptic function must have at 
least two poles (which may merge into a single multiple pole) 
in each of its period parallelograms. From this it follows, in 
particular, that no nonconstant entire function can be elliptic 
(that is, doubly periodic). 

We shall refer to the number of poles of an elliptic function 
in its period parallelogram as the order of the elliptic function. 
The Jacobian elliptic functions sn(t, &), cn(¢t, k), and dn(f, Rk) (in 
particular, the lemniscate functions si¢ and cl?) that we have 
been studying are second-order functions. The results obtained 
by performing any of the four rational operations, namely, 
addition, subtraction, multiplication, and division, on any two 
elliptic functions with the same fundamental periods 2w, and 
20, yields another elliptic function with the same periods. In 
this way, we can obtain functions of arbitrarily high orders. Let 
P(p = 2) denote the order of a function f (4). It can be shown that 
the equation 


fit)=A 


has exactly p roots in any of the period parallelograms of the 
function f, no matter what the value of the complex number 4. 
(Here, multiple roots are counted in accordance with their 
multiplicity.) 

We see that the theorems on the number of roots of the equa- 
tions 


slt=:A, dt=A, sin(f, A) A. cn(t, k)-== A, dn(t, RJA 


in a period parallelogram are special cases of the overall 
principle. 

Without stopping for other properties of elliptic functions, 
let us point out in conclusion that every elliptic function obeys 
an algebraic addition theorem. Weierstrass proved a theorem 
that in a certain sense is the converse of the preceding one: If 
a mermomorphic function 9 (t) possesses an algebraic addition 
theorem and is not a rational function either of t or of an ex- 
ponential function of the form e* (where ais any complex con- 
stant), the function @ is necessarily an elliptic function. (See, 
for example, the book referred to at the beginning of this sec- 
tion.) 
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All the general sines that we have been studying obey addi- 
tion theorems. They are either elliptic functions (as is the case 
with sl¢ and sn(ft, &)) or rational functions of an exponential 
function, as is the case with 


and 


t a ¢ 
z e —-e@ 
sinh t = ————- , 


2 


These facts can be regarded as illustrations of Weierstrass’ 
theorem. 
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